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STUDIES  OF  THE  ELECTROMAGNETIC  CHARACT  E R  1ST  ICS 
OF  MOVING  IONIZED  GASES 

The  Project  Scientist  for  this  task  is  Dr.  Manlio  Abele. 

The  investigations  of  the  microscopic  and  macroscopic 
properties  in  a  weakly  ionized  gas  reported  in  Chapters  II  and  IU 
of  this  summary  report  were  carried  out  under  the  direction  of 
Professors  Piero  Caklirola  and  Luigi  Napolitano,  respectively, 
and  under  tin'  general  direction  of  Prof.  Piero  Caldirola. 

Tice  experimental  investigation  described  in  Chapter  IV  has 
been  carried  out  by  Messrs.  R.  Tomboulian,  M.  Wecker  and 
H.  Medecki  under  the  direct  supervision  of  Dr.  Abele. 

In  Chapter  V,  preliminary  calculations  have  been  conducted 
by  Dr.  Bruce  Gavril  in  the  evaluation  of  the  electromagnetic  pro¬ 
perties  on  the  basis  of  tire  experimental  information  obtained  with 
tie  shock  tube  facility. 
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I 


results  of  the  electron  density  profile  of  the  ionized  argon  behind  the  shock 
for  several  initial  conditions  of  the  driven  section.  In  the  limited  range  of 
frequencies  which  have  been  used  in  these  experiments,  the  results  show  the* 
high  resolving  power  of  the  technique.  In  the  preliminary  part  discussed 
herein,  only  the  case  of  small  absorption  of  the  microwaves  has  been 
analyzed. 

In  connection  with  the  program  of  computing  the  electromagnetic  properties 
of  the  plasma  from  experimental  results,  the  analysis  of  the  propagation  of 
electromagnetic  waves  in  the  nonuniform  medium  inside  the  cylindrical  wave 
guide  is  still  in  progress.  The  analysis  is  carried  out  assuming  that  the 
complex  dielectric  constant  of  the  medium  is  unknown,  and  the  propagation 
equations  are  solved  in  terms  of  measured  quantities  at  the  surface  ot  the 
wave  guide.  The  measured  quantities  correspond  to  the  signals  obtained 
with  the  electric  and  magnetic  probes  which  ate  located  at  the  surface  ot  the 
wave  guide.  Chapter  V  presents  the  results  obtained  assuming  that  the 
electromagnetic  properties  ot  the  plasma  are  constant  in  a  cross-section  ot 
the  tube  and  depend  only  upon  the  axial  distance  from  the  shock.  1  hese 
results  will  be  applied  to  the  calculation  ot  electromagnetic  propel  ties  on 
the  basis  of  the  results  described  in  Chapter  IV.  The  analysis  will  be  extended 
to  the  case  of  nonuniform  transverse  distribution,  assuming  a  suitable  model 
for  the  radial  dependence  based  on  theoretical  calculations  of  the  flow  field 


behind  the  shock. 
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CHAPTER  II 

ELECTROMAGNETIC  PROPERTIES  OF  NON  UNIFORM  PLASMAS 
FROM  MICROSCOPIC  CONSIDERATIONS 

INTRODUCTION 

The  main  object  of  our  research  was  to  examine,  in  a  microscopic  way, 
some  effects  in  a  plasma,  clue  to  inhomogeiicit ies ,  so  as  to  estimate  the 
importance  of  such  effects  on  the  principal  elec  tromagnetic  properties  of 
the  plasma,  as,  for  instance,  conductivity.  In  order  to  reach  such  results 
from  our  view  peint,  it  is  necessary  to  calculate  certain  distribution  functions, 
directly  Reducible  from  Boltzmann's  equations.  If  inhomogeneities  are  present 
(due  to  gradients  of  concentration,  of  tempi  future,  etc.),  as  far  as  we  know, 
no  actually  developed  method  of  calculation  exists. 

For  this  reason  we  decided  to  examine  various  models  beginning  with 
the  simplest  ones: 

j.  Plasma  assimilated  to  a  weak l \  ionized  Lurcntziun  gas,  considering 
electrons  and  neutral  molecules  and  tie  leitinu  Coulomb  interui  tions, 

i.  Three  component  plasma  (elec  trot  s.  ions,  neutral  molecules)  taking 
Coulomb  interactions  into  ai  coui  t. 

Since  the  analysis  of  model  i.  lias  appeared  particularly  difficult  from 
a  general  point  of  view  (for  the  reason  that  it  implies  tne  solution  of  a  system 
of  Boltzmann's  equations),  first  of  all  we  decided  to  carry  out  in  a  detailed 
way,  a  method  of  actual  calculation,  in  order  to  face  the  simplest  model  l. 
even  if  it  is  not  very  realistic.  In  addition,  for  the  scheme  1.  wc  have  decided 
to  examine  separately  the*  various  inhomogeneity  effects  and  also,  to  keep  near 
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to  the  model  studied  from  the  macroscopic  view  point,  we  have  selected  the 
detailed  examination  of  the  influence  of  temperature  gradient  in  a  slab  con¬ 
taining  the  plasma,  whose  walls  are  absorbing  and  emitting  with  a  Maxwellian 
distribution. 

First,  in  Section  L,  the  coipplete  mathematical-physical  characteristics 
(in  our  opinion  not  so  wholly  developed  in  the  literature)  of  modi  l  l  .  are 
described  up  to  the  deduction  of  final  equations. 

Next,  in  Section  <i.,  a  known  analytical  method  used  for  solving  the  equations 
of  Section  l.  (series  expansion  of  eigen  functions  of  t  u  collision  operator)  is 


described. 
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l .  MAT HEMATICAL-PHY SICAL  CHARACTERISTICS  OF  M O DEL  l 

We  consider,  as  in  the  Lorentz  gas  model,  a  binary  mixture  of  a  single 
kind  of  neutral  molecules  and  electrons.  Tii.it  is  to  say,  we  consider  a  weakly 
ionized  gas  (n  «  N)  which  allows  us  to  neglect  (within  tile  limits  of  an  initial 


research  whose  aim  is  to  test  a  niethematic.il  method  for  studying  Boltzmann 
equation  for  a  nonliomogeneous  system)  the  ion  effects,  provided  that  their 
density  is  of  the  order  of  n. 

We  shall  indicate  with  f  (x,  _v,  0  and  F  (x,  V.t)  respectively  the  electronic 
and  molecular  distribution  functions. 

In  the  corresponding  system  of  the  two  Boltzmann  equations  : 

Df  =  .T  (f)  t  J  (f) 
ec  cm' 


DF  .1  (1)  t  J  (F) 

mm  nie' 

we  can  neglect  J  (f)  and  J  (F). 

1  c  me 

Assuming,  then,  the  following  system  for  the  t<  mporal  behavior  of 
f  and  F, 


Df  =  Jcm(f>  n-D 

DF  -  J  (F)  (1.2) 

mm 

we  may  note  that: 

a)  From  the  Equation  (1.2)  the  temporal  behavior  of  the  molecules  is  entirely 
independent  of  the  presence  of  electrons,  while  the  electron  gas  depends  upon 
the  presence  of  molecules,  through  the  collision  term  J 

cm 

b)  The  temporal  behavior  of  the  electron  gas,  given  by  Equation  (1.  1)  is 
affected  by  the  molecular  gas,  as  we  shall  see  later,  only  through  the 
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following  phenomenological  parameters: 

molecular  density,  mean  velocity  and  mean  square  velocity. 

As  a  result  it  is  not  necessary  in  solving  Equation  (1.1)  to  integrate  the 
Boltzmann  Equation  (1.2): 


DF  =  J  (F) 
mtn 


for  we  only  need  the  knowledge  of  the  above  mentioned  phenomenological 


parameters. 

We  shall  begin  by  expanding  the  electronic  distribution  function 


f(x,v,t)  in  spherical  harmonics,  writing: 


f(x.v.t) 


•  r  = '  o" 
‘  o  m  •  e 


°e.n(ji*y*l)  Cem  *  Pem  l)  Se 


em 


(1.3) 


where 


C  =  v°  9  (0)  cos  mO 

em  cm 

S  -  v1’  0  (0)  sin  mO 

em  em 

Q  (0)  being  the  associated  Legendre  functions: 
em 

(m) 

0  (0)  =  sinm0  P  (cos  0) 

em  e 

We  may  write  the  expansion  (1.3)  in  tin-  form 

f(v)  fQ(v)  t  Z_  '  _f,(v)  I  /  (v)  (1.4) 

V 

where  fD  represents  the  isotropic  part  of  the  distribution  function  f,  £|  is  the 

coefficient  of  the  first  anisotropy  and  ^  is  simply  the  remaining  part  of  the 


expansion. 
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Introducing  (1.3)  into  the  Boltzmann  Equation  (1.1)  we  obtain  a 
hierarchy  of  equations,  which  is  easily  obtained  provided  that  we  know  the 
resul'  of  the  application  of  the  linear  collision  operator  J  over  each  term  of 
the  expansion  (1.3). 

To  this  end,  we  may  note  that  the  operator  J,  in  an  approximation  for 
which  one  has  simply  77  =  0  (perfect  Lorentz  gas),  possesses  the  following 

ivl 

properties,  for  each  isotropic  function  n(v): 

J_(a)  =  0 


Jem(aCem>  --  -  VCcm 

.1  (aS  )  t  -  V  «S 
cm  em  e  cm 

In  the  imperfect  Lorentz  gas  scheme,  we  shall  have: 


J  (n)  =  O,  (  — ) 
em  M 


(1.5) 


J„ni{aC^J  ^  -V  ,iC  1  ♦  Oz  (  — ) 

em'  cm'  em  _  ‘  M 


J,.M1(aS  )  -  t>S  [l  ♦  0,  ( 

cm'  c  m  e  em  l_  5 


M 

— ) 
M 


(1.6) 


where  we  have  used  the  Landau  symbol  Oi  (x),  meaning  a  quantity  which  is 
simply  of  the  order  of  x  . 

The  quantity  O  (  — )  has  been  calculated  by  others  (  2)  under  the 

1  M 

hypothesis  of  a  Maxwellian  molecular  distribution  F,  obtaining: 


.  ,  .  _  m  1  3 

J  (fl)  =  -  —r-  - — 

cm'  M  V2  3v 


,,  3  ,  kT  9a. 
P  v3  ( a+ - ) 

mv  ov 


(1.7) 


On  the  other  hand,  we  have  proceeded  to  the  evaluation  of  the  quantity  Oj(m/M) 
keeping  an  entirely  arbitrary  (and  hence  generally  anisotropic)  molecular 


distribution. 
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The  results  of  our  calculation  is  given  by: 

i  a 


J  (a)  = 

em  v 


.  J 

2  3  v 


A(v)  4  B(v)  •  V  +  C(v) 

V  —  2 

v‘ 


:  VV 


J 


where: 


A(V)  -  H  v3  Q  +  1/  v5  3 


V2  3q 


3  v 


(1.8) 


B(v)  l/{  v2a 

CM-  (.  <  3 


i  i  a  v 

IV  -  -  V  )  v  -  -  — i  v2  |  a 

1  2  2  6  3  v 


(1.9) 


We  have  followed  Chapman  and  Cowling's  notation,  in  writing 

_  3 


V  V 


^  (v. v  -  —  v2  6  )  V.V. 


1  j  3 


ij  J  i 


IJ-Tl 


while  and  t2  arc  the  first  two  r<  luxation  frequencies  1 


1<J (\ )  -  2"  Nv  ,  s in  0 


1  -  P,  (cos  0) 


u(0,  v)  dO 


1/  (v)  =  2 72  Nv  sin  0 


1  -  P  (cos  0) 
2 


(7(0,  v)  do 


(1.10) 


C(0)  |v  -  V|)  being  the  differential  cross  section  for  electron-molecule 

collision,  in  a  frame  of  reference  fixed  with  one  of  the  two  particles.  The 

determination  of  .1  (nC„,„),  J  (oS  )  has  been  done  disregarding  the 

cm  em  cm  en* 

02,  03  terms. 


A  general  account  of  the  method  followed  in  the  derivation  of  formula 


(1.8)  is  given  in  Appendix  1. 
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Taking  into  accout  that  £  is  an  isotropic  function  and  is 

O  v 

simply  a  linear  combination  with  isotropic  coefficients  of  C,0  ,  Cu,  S(l  , 

v 

we  are  now  able  to  calculate  Jem  (f0).  Jcm  ^ 

If  we  perform  a  further  approximation^)  forgetting  the  tcrm^(v) 
in  (1.4),  we  at  last  arrive  at  the  following: 


Df  = 


J 


em 


(f) 


_L  — 

v2  Sv 


r 

V  V  V 

'  A(v)  t  n  (v)  -=—  *  V  +  C(v)  =^=- 

v2 


I 


J 


V 

i 


v 


V 


(1. 1 1) 


where  A(v),  B(v)  and  C(v)  are  given  by  (1 .9). 

The  calculation  of  Df  is  immediate,  finally,  by  means  of  the 
orthogonality  properties  of  the  spherical  harmonics  functions  on  the  unit 


sphere,  we  obtain  the  following  system: 


-S 


dt 


at 


c 

_b 

1  ^  ■ 

3  m  v 2 

d\ 

M 

V2 

Sfo 

3m  v 

Sv 

e 

3fo  E 

;  t  - 

m 

"V  “ 

tS  f  , 

o 

,  f  . 

1 

2  ,  ,  >  m  1 

r  h. f  i  )  =  —  —7  — 

M  v2  5  v 


v3(fo 


P.  12) 


1  * 


♦  V?fo*  -  *  *  mc<iL^)=  ^  rv  ^  fu> 


)  As  usual  in  this  kind  of  reasoning,  we  may  evaluate  "a  posteriori1  the 
usefulness  of  this  approximation,  by  calculating  explicitly  the  difference 
between  the  results  obtained  and  those  based  on  the  successive  approxi¬ 
mation,  or  "a  priori1'  we  may  introduce  some  plausible  arguments,  like 
those  used  by  Davydov  (  1  )  and  Gurevitch  (  4  )  . 


APPENDIX  1  -  DERIVATION  OF  EQUATIONS  1.8  AND  1.9 
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Let  n s  state  here  some  formulas  which  we  shall  use  later  on, 
directly  declucible  from  the  dynamics  of  elastic  (classical)  collision  between 
two  particles  (in  our  case,  electron  and  molecule). 

Let  us  consider  an  elastic  elec  iron-molecule  encount*  r ,  with 
velocity  v,  V  respectively,  well  before  the  collision. 

Let  us  put: 

j.  -  V  -  v 

where  g  is  the  relative  velocity  of  approach  of  the  molecule  in  a  frame  of 
reference  fixed  with  the  electron  (see  Figures  1  and  Z  at  the  end  of  this 
Appendix). 

0  is  the  scattering  angle;  from  Figure  Z  wc  hive  at  once: 

0 

g-g'*.£gsin  —  k  =  Z  ( it  •  )  k 

where  k  (apse  line  unit  vector)  represents  tin-  clirecli  >n  of  the  external 
bisectrix  of  the  angle  between  the  asymptotes  of  the  relative  molecular 


trajectory. 


We  have  thus: 

y 

V  -  v  =  V'-v'  *  Zg  sin  —  k 

from  which,  by  momentum-conservation: 

m 

V'-V  -  -  —  (v*  -  v) 

Z 

v'  -  v  -  -  (g  •  k  )  k 

-  1  ‘  A  -  -  - 

ZA  ,  ,  .  , 

V'  -  V  =  -  -  (g  •  k)  k 

—  —  1  1  A  - - 


(A  1  -  1 ) 


1 1 


having  stated  A  =  — 

M 

In  the  frame  of  reference  sketched  in  Figure  3,  we  see  that  the 

components  of  k  are 

,  0 

k  =  -  cos  —  cos  O 

*  2 

,  0 

k  =  -  cos  —  smC 

y'  2 

k  =  sin  - 
2 

We  pass  from  this  frame  to  the  one  sketched  in  Figure  4  through: 
n  -x'  sin  (3  +  y'  sin  (  cos  (3  *  •/.'  cos  f  cos  |3 
y  -  x'  cos  (3  +  y'  sin  f  sin  |3  t  /. 1  cos  c  sin  (3 
■/.  -  y'  cos  f  -  sin  C 
when  (see  Figure  1): 


v  V 

sin  (  -  —  -  —  cos  u 

K  R 

V 

cos  f  -  sin  n 

R 

So  we  can  calculate  the  components  if  k  in  the  frame  of  r<  ference  of 
Figure  4.  Using  the  first  of  (A  1-1)  we  have  the  following: 


v'  (1  f  A)  =  - 
x 


(v  -  V  cos  n)  sin  U  sin  O  -  V  (1-cos  0)  sin  ci 


cos  (3 


)  g  sin  0  sin  (3  cos  <3 


v'  (  1  +  A)  = 

y 


-  (v  -  V  cos  a)  sin  U  sin  O  -  V  (1-  cos  U)  sin  a 


sin  (3 


(A  1  -  2) 


g  sin  0  cos  |3  cos  O 


v'  (14-A)  =  v  (A  t  cos  0)  i  V 

z 


(1-cos  0)  cos  a  -  sin  Q  sin  a  sin  C 


So,  squaring  and  summing: 


12 


v'  = 


1  4  A 


'  1  4  2A  cos  0  i  A2  4  2B2  (1-cos  0) 


L 


1  2B 

where  we  have  put: 


( A  —  1 )  cos  q(  1  -  cos  0  )  -  vAt  1)  sin  a  sin  0  sin 


V 

B  =  — 
v 


We  now  treat  A  and  B  as  infinitesimal  quantities  and  eliminate 
2  3  3  4 

quantities  which  contain  A  ,  A  ....  B  ,  B  . AB, 


(A  1-  3)  becomes  : 

1 1  A  cos  Q  -  A  +  B2  (1-cos  0)  -  B 


V'  -  V 


cos  a  ( 1 -cos  0) 


4  sin  0  sin  a  sin  W 
To  the  same  order: 


B2 


cos  a  (1-cos  0)  t  sin  o  sin  0  sin  C? 


v'  -  \  v  (B2  -  A)  (1-cos  0)  -  B  cos  a  (1-cosO) 

i.  L 

B2  1  ^ 

i  sin  a  sin  Q  sin  O  -  —  l  cos  n(  1  -cosO)  1  sin  u  sin  0  sin  O 


(v'-v 


r  .  i2 

-  v2B2|^cos  a  (1-cos  0)  t  sin  a  sin  0  sin  O 


m  ,  b2  .  2 

=  v  j  1  *  —  sm  n  -  B  cos  a 


g  -  v  1  1  B2  -  2B  cos  a 

L  j  u  *■ 

|  B2  da  B2  ,  , 

a  (0,g)  -  a(0,  v)  t  v  j  —  sin2a-  Bcosa  j—  4  —  v  cos2a 


0  (0,  g)  g  =  a  (0,  v)  v-Bcosa  (v  - —  t  a )  4  O  (B2) 

dv 


We  must  now  simply  evaluate  the  behavior  of  the  operator  Jom  over  an 
arbitrary  isotropic  function  q(v),  that  is  the  quantity: 


(A  1  -  3) 
1/2 


(A  1-1) 

(A  1-3) 

(A  1-6) 

(A  1-7) 

(A  1-8) 

(A  1  -9) 


JemM  =  \(a'F'  -  Q  F)  Sa(°-g) 


(A  1  -  1 0) 


1  5 


The  proceeding  collision  term  has  been  evaluated  by  noting  that  in  each 
electron-molecule  collision,  the  variation  of  the  modulus  v  of  the  electron 

i 

velocity  is  given  by:  - 

v'  -  v  -  O  (  —  )lU  (A  1  -  1 1 ) 

M 

(see  (A  1-5)  ) 

As  a  consequence,  the  integral  operator  J  can  be  put  in  a 
differential  form,  following  the  same  technique  used  in  the  transformation  of 
the  Boltzmann  elastic  collision  operator  into  the  Fokkcr-Planck  form 
(references  (  5  )  (6  ) ) 

Let  us  introduce  a  suitable,  but  largely  arbitrary  isotropic  function 
t  (\);  according  to  a  general  well  known  r  ile:^  ^ 


J  ^  (n)  /  (\ )  dv  =  |^(v')  -  </(v)  j  g  n(v)  F0(0,  g)  dfldV  dv  (Al-12) 

Furthermore,  retaining  only  terms  of  the  order  of  m/M,  we  have  according 
to  (A  1  -  1 1 ): 

^  i  _£  V 

(A  1  -  1  3) 


(%  *)  -  )  (v'-\)—  t  i  (v'-v)2  ~~ 

2  ev2 


Inserting  ( A  1  -  1  3 )  into  (Al-12)  and  integrating  by  parts,  we  find: 


1  o 

J  (a)  -  -j- 
em  v‘  ev 


-v2  M(v)  ft  (v)  t  —  — (v2X(v)a(v)) 
1  3v 


where 


M(v)  =  (  0(0,  g)  gF  (v'-v)  tin  dV 
X(v)  0(0,  B)  gF  (v'-v)2  dOdV 


( A  1  -  1  -4 ) 


(A  1  -  1 5) 


(*)  See,  for  instance  references  (2)  ,  (7  )  . 
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So  f;ir  as  we  choose  $  (v)  to  satisfy  the  following  relations: 


|nM^  v2 

a) 

\  LL  2 

aA  .  \c 

0  V 

a)  1  * 

v-o 

-  J 

V  O  U  V 

In  a  way  consistent  with  the  approximations  used,  the  two  preceding 
functions  U(v)  and  X(v)  have  to  he  evaluated  to  the  first  order  in  the  quantity 
m  /  M . 

To  the  first  order,  C(Q,  g)g(v'- v)  is  given  by  formulae  (Al-S)  and 
(A  1  -9)  as  follows: 

0g(\  ’ - \  )  =  0(0,  v)  v2B2  ( 1-cos  0)  -  0(0,  \  )  v2A  (1  -cos  0) 


-  CT  v2  B  ios  a  (1-cos  0)  -  0v2 


B * 
2 


cos2u  ( 1  -cos  0)2 


-  °v2  —  sin20  sin2n  sin2C  *  B2  cos2a  v2  (v  —  i  0)  (1-cos  0) 

ov 

where  quantities  containing  sin  O  have  been  overlooked,  because  a  successive 

integration  is  performed  over  O,  from  o  to  217  .  From  the  definition  (A  1-1  5)  of 

we  (  an  obtain  it  from  the  preceding  relation,  byr  integration  over  0  and  O 

and  over  all  the  velocity  space  V. 

Leaving  aside  the  (direct)  integration,  we  have: 


M  (v)  *  —  V 
v 


A  V  i  v  -  v  V 


i 

2  v 


b 

2 


cU' 

I 

J  5  V 


(A  1-16) 


where 


77 

=  277  Nv  j  sin  0  (l-cos  0)  0(0,  v)  dO 

o 

77 

a  =  277  Nv  sin  0  (1-cos  Q)2  0  (0,  v)  dO 
o 

TT 

b  =  ZTT  Nv  ^  sin3  0  (J  (0,  v)  dO 


o 
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and  having  stated: 

V  =  —  \  dV  V  cos  a  F 
l!  N  )  - 

V2  =  —  '  dV  V2  cos2n  F 

II  N  )  ~ 

V2  =  —  dV  V2  sin2n  sin^  F 

i  n)  - 

In  the  same  way,  to  the  first  order  in  m/M,  CJ  (J,g)  g  (V  -  ' )  following 
formulae  (Al-f  )  (A  1-9)  id  given  by: 


CT(O-g)  g(v'-\)2==  dv  V2cos2n  (1-cos  0)2  I  O  \  V 2  sin2  0  sin  n  sin^ 

overlooking  quantities  containing  sin  O. 

Following  the  definition  (A  1-1  5)  we  ha  e  therefore: 


X  (v)  = 


* v* 


( A 1  -  1 H ) 


Let  us  now  introduce  the  quantity: 


V  -  2?TNv 
z 


77  1 

\  sin  l)  1 

cos  20  -  —  ) 

2 

o 

L 

J 

o  (0,  v)  do 


(A  1-19) 


(A  1-18)  can  he  written: 


A  (v)  -  UVt  -  VJ  V2(  .  ±  vt  V2  , 


(A  1  -  20) 


with 


V7  t  V2  =  V2 

II  l 

Introducing  now  (A  1-16)  and  (Al-20)  into  (A  1  - 14),  and  leaving  aside  the 
direct  calculations,  we  have: 


J  (n) 

em 


V2 


_3_ 

d  v 


{ 


A 1  (v)  +  IV  (v) 


1  C’(v) 


I 
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where: 


=  —  t/,  v!atv(il'  -  )  V2  a  +  j*  v2  V2a  +  1  v*  ^  ^  V2 

2  2  6  6 


A(  v ,  - 

M 


"  B'  (v)  l/,  v2  a 


i  1  ,  *vz 

C'(v)  -3v  (i  V  -  Vt)  a  -  ~  v*  —  u  -  v 


1  5a 

1  <  -  V  -  l/j)  — 

2  2  5v 


Furthermore,  by  noting  that: 


Vj,  =  f  •  V 


V 2  -  i  V  v  :  V  V 


and  puting: 


we  have  at  last: 


C  1  ?  X 

V  V  :  V  V  -  ~  V‘  O 


1  5 


v®v  _ _ 

T  (n)  ,  -i-  -^-4A(v)  .  »(v)  4-  V  «C(v)=-:  V  V  ’  (M-21) 

'Mil  '  ’  v2  ?v  v  —  v2 


where 


~\ 


m  i  j  v  a  a 

AW=  -  *,»■••*»*  7  7 


E3(\)  -  V.  v2  n 


(A  1-22) 


■  5  a  1  1  i 

C(v)  [u^  -  -  V)  v2  ^  <  3  L(«/,  -  -  ^  )  v  -  -  v 


>  2  2  ov  L 

that  is  formulae  (1,8),  (1.9)  of  Section  1, 


(A  1-22)  can  be  put  in  the  form: 


1  5A_  1  5B  vv 

Jem(a)  =  ^  3v  +  ^  5v  v* - * 


(A  1  -  23) 
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At  this  point,  we  are  able  to  calculate  the  quantity: 


Jem  <f> 


where  f  is  the  electronic  distribution  function: 


f  =  f0(v)  (v) 

v  — 

Taking  account  that  Jem  is  a  linear  operator,  we  have: 


J  <f)  -  J  (fG)  -  V-  •  f. 
em  em  °  1  v  J 

by  using  the  Bayet-Delcroix-Denis.se  approximation  (see  Reference  (1)  )  , 


One  can  evaluate  )  from  (Al-23)  by  simply  substituting  n  with  f  and 

a r r i \ e  .it  the  formula  (1.11)  of  Section  1. 


18 


FIGURE  I 


.  plane  through  the  origin,  perpendicular  to  the  inconiing 
'=  half-plane  terminated  by  the  straight  line  n,  ovir  v. 

molecule. 

■  "=  reference  half-plane 


b  =  inpact  parameters 
0=  a/.imuthal  angle 
0  =  scattering  angle 
g1  =  relative  molecular  velocity 


after  encounter 


incoming 


(This  Figure  is  intended  to  be  traced  on  the 
T1  half-plane  of  the  preceding  Figure  1;  k  is 
the  unit  vector  of  the  so  called  ’apse”  line.) 


\ 

\ 


21 


ANALYTICAL  METHOD  FOR  SOLVING  EQUATIONS 

We  now  approach  the  study  of  the  system  (1.  12)  limiting  ouselves  to 


the  stationary  state  when  there  is  no  magnetic  field. 

It  is  thus  possible  to  express  f,  (x,  v)  as  a  function  of  f^  (x,  v)  and  to 
reduce  the  system  to  only  one  differential  equation  in  fQ  (>b  v ) 


I  ? 
3  _ 


_v  o  f  . - 

--  l‘o  mV 


of  1 

_ o  E  t  — jjj 

5  v 


?  — 

—  (y,  v2  f  )  V  ✓ 
dv  v  1  o'  - 

j 


(2.  1) 


m 


3mv2  3v 


-  ^  E  •  £.  >  =  M 


& 

b  v 


y,  vJ  (fn  » 


MV2  »f0\ 

3,nv  77'j 


Then  we  studied  this  elliptic  partial  differential  equation  in  the  variables 
x  and  v. 

As  we  said  before,  the  presence  of  the  molecular  gas  appears  only 
through  the  mean  molecular  velocity  J  ,  Vz  and  the  mol,  cnlnr  concentration  N. 


We  developed  f()(x,v)  in  a  series  expansion  of  eigenfunctions  of  the 

&  f  /  MV  2  o  f 

3  if  ,  -  - o  In  doing  so  it  was 

•  Vl«  3m  V  5  v  j 

necessary  to  specify  the  type  of  electron-molecule  interaction  we  were 


m  j _  _ 

operator  .1  -  v2  }v 


conside  ring . 

We  have  considered  two  cases  corresponding  to  central  interaction  of 
the  Maxwellian  type  and  the  rigid  spheres  type  (see  Appendix  2). 

In  the  first  case  the  eigenfunctions  of  J  are: 


J 


mv‘ 


mv‘ 

2k~T 


(1/2) 


111V 

2kT 


l,-j  -c  d/2) 


i  e  L 


(0 

J 


(2.  2) 


(where  £  —  —  is  the  non-dimensional  electron  kinetic  energy);  in  the 

2kT 


second  case  they  are: 
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mv* 


<.  2  kT 


<n  f 


mV 


N 

)  J 

r 

’  1 

(2.3) 


Thus  wc  have: 

a)  For  Maxwellian  interaction: 

oo  (1/2) 

f  =  .  a  (x)  e  L  (f) 

o  X  11  -  n 

*n  =  0 

b)  For  rigid  spheres  interaction: 

oo 


(2.4) 


f  - 

o 


n  ( x )  c 
n  — 


-f 


(1) 


L  (f) 
n 


(2.5) 


V-TT 


(k) 


Using  the  orthogonality  properties  of  the  L(1  polynomials,  we  arrive  at 
two  systems  of  infinite  differential  equations  in  an(x)  and  nn(x)  which,  if 
the  electric  field  is  neglected  and  V  0,  takes  the  form: 

case  of  elastic  -  spheres  interactions:  (2.ti) 

of  (an(x)  )  '  (n  t  2)  ^(nn(x))  r  (n  •  2)  (n  t  3)C  (nn(x)  )  - 


-  n$(«n_j  (x)  )  -  2n  (n  •  2)£  (n^  _  ,(x)  )  ♦  n(n-  1)£  (an_  ^(x)  )  = 

=  3  V2  Nz  ~  nn  (x) 
o  M  n 


(n  0,  1,  Z,  .  . ) 


1 3 


Case  of  Maxwellian  interaction: 


(2.7) 


!-n(n.l)(n-2)C(an.3(x)  )+  n(n-  1  >j8  +  n(n"  1  *  (4n*2  ^  (an- 2(-*  1  ' 

-  .  9 

'!  ni(nn-l)  *  n<3nf2)^  (nn-l)  ‘  n(n  f  |  )  <6n  '  2*  C  (an.j)  1 


l 


{ 2n  +  |  )£  (an)  ♦  (n  4  \)  (3n  »  |  )£  (an)  +  (n  t  |)  (n  t  |)  (4n  J)  £  (■'„) 

(n  +  -)<£  (a  , )  t  (n  4  2)  (n  *  -)^(a  )  *  (n+  |)  (n  f  -  )  (n  *  -)£ (a  ,  j) 

2  iH  1  2  2  nt  1 

N2 


J 


n  an(x) 


( n  —  0.  lf2»...) 


MkT 


where  : 


f 


L 


oy  s  72  -  V  log  N  ’  V 


^  ==  2  V  log  T*  7  *  V  2  lot-  T  -  (  V  log  l')2  -  V  log  N*  7  log  T 


£  =  (V  log  J  )2 


In  the  Maxwellian  interaction  ase  we  have  also  evaluated  the  contri¬ 
bution  duo ,  s<  pa rate  1>  ,  to  the  presence  of  both  an  electric  field  and  of  a  mean 


molecular  velocity. 
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APPENDIX  2  -  EIGENFUNCTIONS  AND  EIGENVALUES  OF  THE 
COLLISION  OPERATOR  J 
a)  Maxwellian  interaction  (*) 


6mV2  1/2  <1  ,  3/2  ,  do 

- —  f  77  f  (Ct  elf  >1 


M' 


M 


We  consider  the  following  eigenvalue  problem: 


(A  2-1) 


or 


JM(C  )--  -  \o 


Zf  02.  M2C  ■  3)  -  •  m-ML-  io  =  o 
df2  d€ 


(A2-2) 


(A  2- 3) 


Introducing  the  new  function 


+  < 

(()  e  O  (c) 


(A 2- 3)  becomes: 


cl 2  b  }  ,  d .  MX 

f  —  *  (  -  -  o  —  •  - - f 

2  dC  6mr 


(A2-4) 


2  dC 

o 

Imposing  on  *(<).  as  boundary  conditions  in  (-space,  analyticity  m  f*  0 
and  polynomial  behavior  for  f  >m,  the  solutions  of  (A2-U  are  the  Laguerre’s 
generalized  polynomials  of  order  1/2: 

,<w>lo.  •"  >  | nl  1/2 ;  - 

Ln  ^  ,  n-U  V  I 

'v  =0 


The  maxwellian  electron-molecule  interaction  model  corresponds  to  a 
repulsive  force  varying  as  the  inverse  fifth  power  of  their  distance  apart. 
In  this  case  the  collision  frequency  does  not  depend  on  the  electron 
velocity:  V  -  VQN. 
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corresponding  to  the  eigenvalues 


6ml'  2 

_ o 

M 


(n  =  0,  1,2...) 


(A2-5) 


b)  Rigid  sphere  interaction  (*) 


6ml'  _1  d  dp 

J  ‘  *  |'  (°*  “> 
K  M 


We  consider  the  following  eigenvalue  problem 

J  (O  )  -  -  Ac 
R 


or 


MX 


<2  ^SL  f  (2tO  —  +  (2  t  ■  7  )  o  -  0 


d«* 


df 


Introducing  the  new  function 


(A  2- 6) 


(A  2-8) 


MO  c  O (  C  ) 
(A2-8)  becomes: 

d2t 


e  -  *  (2- f )  —  t 

d<2  df  Gn^2 

o 


*  =  0 


(A2-9) 


Imposing  on  .  (»)  the  same  boundary  conditions  in  f- space  as  in  the  Maxwellian 
case,  the  solutions  of  (A2-9)  are  the  Lnguerre's  generalized  polynomials  of 

order  1,  I 

I  n 


(f)  corresponding  to  the  eigenvalues 
(n  0,1,2,..) 


6m  y  2 

*n  =  _ 1  ” 


(A2-10) 


M 


(*)  In  this  model  we  suppose  electrons  and  molecules  are  rigid  spheres.  In 

this  case  the  collision  frequency  depends  linearly  on  the  electron  velocity 

V-  V  N  v  . 
o 


1 

Zb 

I 

SYMBOLS 

I 

light  velocity  in  vacuum 

k 

■ 

Boltzmann  constant 

■  m 

electron  mass 

1  c 

electron  charge 

r 

1 

electron  (classical)  radius 

*  n 

electron  density 

1  - 

electron  velocity  (modulus  v) 

M 

| 

molecular  mass 

■  R 

molecular  radius 

I  - 

molecular  velocity  (modulus  V) 

1 

N 

molecular  concentration 

J 

T 

molecular  Temperature 

1 

V 

mean  molecular  velocity 

I 

V2 

mean  square  molecular  velocity 

! 

E 

electric  field 

j 

H 

magnetic  field 

D 

differential  Boltzmann  operator 

d 

D—  —  tv  •  V  t  x  ’V  v 

-  at  -  -  -  -- 

i 

1  J 

integral  Boltzmann  operator 

V-  i^Nv  =  TT(R^r)2  Nv  collision  frequency  for  rigid  spheres  interactions 

V^V'  N 
o 

^collision  frequency  for  Maxwellian  interaction 

i 

1  A=  — 

*  dx 

constant  temperature  gradient. 

1 

i 

knudsen  number 
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CHAPTER  I_H 

STUDY  OF  ELECT ROMAGNETIC  PKOPERT 1 ES  OF  NONUNIFORM 
PLASMAS  IN  THERMAL  EQUILIBRIUM,  ISOTROPIC  CASE 


INTRODUCE  ION 

In  this  chapter  results  ol  the  study  of  the  electromagnetic  properties  of 
a  nonuniform,  macroscopica  Ily  neutral  plasma  ar  reported. 

The  line  of  approach  is  a  macroscopic  one,  based  on  the  trie rmodynamic s 
of  irreversible  processes. 

A  three  fluid  model  i.*,  assumed;  tnc  plasma  is  considered  to  be  at  rest 
in  a  suitable  reference  ir.iim,  t  ie  components  an  in  mutual  thermal  equilib¬ 
rium  and  imposed  niacin  tie  fields  are  absent.  I'he  nonuniformities  are  those 
connected  with  the  presem  e  of  gradien  ts  of  stale  parameters  such  as  tempera¬ 
ture,  density  and  i  onm  t  rations . 

The  presence  of  t  u  se  gradients  indu  es  an  elec  trie  field  (i<  e.  f.)  in  t  ie 
medium,  even  in  t  a  absence  of  external  ■  le  tri  liedds.  It  will  je  shown  that, 
in  the  stationary  state  and  in  condition  ol  mac  roscopi  neutrality,  the  induced 
i  .  f.  depends  only  on  one  independent  gradient  \v  i*  ii,  for  convenience,  is  taken 
to  tie  that  of  temperature.  An  expressi  m  for  the  i.e.f.  is  derived,  which  is 
valid  for  a  most  gcncr.il  medium  and  which  is  later  simplified  to  the  ease  of 
an  imperfect  Lorentz  gas  (i.e.,  weakly  ionized  plasma  for  which  the  mass  r.itio 
between  the  negative  charge  and  neutral  mol  ill  s  is  much  smaller  than  one). 

The  electrical  conductivity  0  defined  as  the  ratio  between  the  electric 
current  I  and  the  electric  potential  gradient  in  a  first-order  stationary  state 


in  which  a  given  nonuniform  distribution  of  V  §  is  maintained  on  the  boundaries 
of  a  system  is  also  determined. 

The  phenomenological  coefficients  are  first  expressed  in  terms  of  thermo¬ 
dynamic  properties  of  the  plasma.  It  is  shown  how  these  coefficients  can  be 
expressed  in  terms  of  only  simple  binary  molecular  diffusion  coefficients  and 
thermal  diffusion  coefficient  s  pertinent  to  the  plasma  constituents. 

This  last  step  makes  it  possible,  by  utilizing  the  available  experimental 
and/or  theoretical  data  on  electron-ion,  elet  t  ron-molecule,  and  ion-inolecule 
collisions,  to  express  the  plienomcnologi  .il  ■  oetiicients  in  terms  of  the  local 

thermodynamic  state  of  the  medium. 

These  results  provide  the  necessary  pri  lii  inary  information  needed  to 
proceed  to  the  actual  solutions  of  fluid  dynami  .ml/or  electromagnetic  phenomena 
in  a  nonunifo rtn  plasma.  In  addition,  they  already  lend  themselves  to  interesting 
order  of  magnitude  analysis  on  the  el'fe<  1 .1  indu  ed  by  the  nonuniformities,  The^ 
can  be  used,  for  instance,  to  define  rai  g<  wit  ii  which  these  inhom ogeneties 
can  be  neglected  and  t  e  medium  treated  as  homogeneous,  and  to  lend  theore¬ 
tic. il  substantiation  to  experimental  findings  During  tie  last  two  years,  for 
instance,  some  experiments  were  made  by  S.  Klein  U  Sac  lay  (France)  on  the 
direct  conversion  of  thermal  energy  into  electrical  energy  (Reference  15). 

The  experiments  were  performed  with  a  glass  bill  in  which  two  metallic 
electrodes  were  immersed  in  an  ionized  mercury  gas  at  rest.  1  he  electiodes 
were  kept  at  two  different  temperatures  and  the  ionization  was  obtained  through 
an  H.  F.  source.  As  a  consequence  of  the  temperature  gradient  an  induced 
electric  field  was  detected  and  measured  as  a  potential  drop  between  the  two 
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electrodes.  It  was  possible  to  observe  th.it  the  i.e.f.  increased  with  the 
temperature  difference  and  the  degree  of  ionization.  These  results  qualita¬ 
tively  agree  wit.,  the  theoretical  conclusions  wh.ich  are  arrived  at  in  the  present 
report. 

The  steps  of  the  analysis  herein  reported  are  as  follows:  In  Section  l., 

Part  l.,  the  irreversible  thermodynamic  description  of  the  system  is  performed. 

I  lie  extensive  and  intensive  state  parameters  and  the  pertinent  mass  and  energy 
fluxes  for  a  mixture  of  three  thuds,  oiu  with  negative  and  one  with  positive  charges, 
art  defined.  1  hi  dynamic  relations  between  fluxes  and  generalized  forces  are 
established  and,  throng. i  suitable  use  ol  the  bn  si.  theorems  of  the  thermodynamics 
<i  irreversible  proi  esses,  the  general  expression  for  the  mass  fluxes  in  terms 
o:  electron  concentration,  temperature,  and  .  lei  trie  potential  gradients  is 
arrived  at. 

In  Section  I.,  Parts  Z.  and  -1.  ,  the  general  expressions  for  tne  i.e.f.  and 
the.  electrical  conductivity  are  obtained. 

In  Section  I.  ,  Part  the  problem  ol  r  lating  these  phenomenological 
coefficients  to  binary  transport  coefficients  is  considered.  By  suitable  trans¬ 
formation  ol  fluxes  and  aftinitii  s  ii  is  shown  how  they  can  be  expressed  in  terms 

of  three  binary  diffusion  coefficients  Du,  IB,,  l),,  and  two  thermal  diffusion 
T  T 

coefficients  D,  ,  Dg  ,  which  refer  to  the  plasma  constituents. 

In  Section  A,  Part  l.,  the  evaluation  of  the  molecular  and  thermal  diffusion 
coefficients  foi  a  plasma  is  carried  out  and  in  Part  Z.  an  analysis  of  the  order 
ol  magnitude  for  the  transport  coefficients  is  performed  for  the  case  of  a  weakly 
ionized  imperfect  Lorentz's  gas. 
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Section  3.  deals  with  an  actual  determination  of  the  order  of  magnitude 
of  the  phenomenological  coefficients,  with  discussions  of  their  properties  and 
with  their  actual  evaluation  for  indicative  cases. 


SECTION  1 


l.  1.  IRREVERSIBLE  THERMODYNAMIC  DESCRIPTION  OF  THE  MEDIUM 

We  consider  a  plasma  formed  by  electrons  (subscript  l)  ions  (3)and  neutral 
molecules  (3)  in  mutual  thermal  equilibrium  (i.e.  ,  T]  =  T>  =Tj=T).  It  is 
further  assumed  that  the  plasma  is  at  rest  in  a  suitable  reference  frame  and 
that  no  imposed  magnetic  field  is  present. 

It  is  of  basic  importance,  for  a  clear  understanding  of  the  physics  of  the 
phenomena  being  studied,  to  proceed  to  a  rigorous  irreversible  thermodynamic 
description  of  the  system.  This  procedure  will  also  help,  on  one  hand,  to  point 
out  the  number  anti  types  of  simplifying  assumptions  which  will  be  made  and, 
on  the  other  hand,  to  have  a  clear  overall  picture  of  the  kind  of  results  that 
can  be  obtained  with  this  approach. 

The  fundamental  extensive  thermodynamics  variables  for  the  subject 
system  are:  the  intt  rnal  energy  U,  the  volume  V,  the  masses  M,,  M,  ,  and  M, 
of  the  three  species  and  the  total  electric  positive  and  negative  charges  E(  and 
E,.  present  in  the  volume  V.  The  fundamental  relation  S  =S  (V  )  (see  Reference. 
33)  cun  then  be  written  as; 

3  L 

TS  =  U  +  pV  Hi  Mj  -  '  f.  E.  (1) 

i=l  T=l 

where  S  is  the  total  entropy;  T  the  temperature-,  p  the  pressure  of  the  mixture. 
The  quantities  gj  and  fj  are,  respectively,  energies  per  unit  mass  and  energies 
per  unit  electric  charge  and  are  thermodynamically  defined  in  terms  of  the 
partial  derivatives  of  S  with  respect  to  the  corresponding  extensive  variable. 
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We  make  the  following  fundamental  hypotheses: 
Ei  =  ci  Mi 

Bj  (XR)  =  Pi  (U  V.  Mj) 
fi  (Xk)  =  0 


That  is: 

The  total  electric  charges  are  proportional  to  the  total  masses  M,  and  M.  . 
This  implies,  among  other  things,  t  lat  then-  is  only  one  type  of  negative  and 
positive  charge  carrier. 

The  quantities  g;  coincide  with  the  chemical  potenti.il  pj  pertinent  to  the 
mixture  considered  in  the  absence  of  electric  charges  (cq  =0), 

The  quantities  f  are  independent  of  the  extensive  variables  Xj  ami  thus 
coincide  with  the  electric  potential  6  . 

It  can  he  shown  that  these  two  latter  hypotheses  are  equivalent  to  the 
hypotnesis  of  separability  of  the  energ\*  U,  in  the  sense  that  the  total  energy 
of  the  system  is  simply  the  sum  of  the  "mechanical  energy"  of  the  noncharged 
system  for  given  S,  V,  M-  (e-  =0)  plus  the-  V  It  iric.tl  energy"  (cq  \1,  t  cw  M,)0. 
This  hypothesis  is  certainly  valid  in  the  subject  case  (cfr.  Reference  19). 


Equation  (l)  can  then  be  written  as: 

s  “  T*  fv'Z 


U) 


i=  l 


where  |ij  is  the  electromechanical  potential  defined  as: 

Pj  =  P,  +  "i  «t> 

and,  obviously,  p  ,  =  p4. 


(3) 


*  A  rigorous  statement  would  involve  the  concept  of  free  energies  (cfr.  Ref  19). 
It  can  be  shown,  however,  that  if  the  "free  energy"  is  separable,  so  is  the 
internal  energy. 
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If  and  Ji  (i  *  l,  L,  3)  are  the  energy  and  mass  fluxes,  respectively, 
using  the  basic  assumption  of  the  thermodynamics  of  irreversible  processes 
one  writes: 


ju  =Loo*  f  M-Tr1-'  ■¥  +  uo>*l-X£- 


j.  =l-.oV  (y>  +  Ln  V(  Y-}  L> H) 


'I  lie  last  relation  follows  from  the  mass  conservation.  The  quantities  L  — 
are  the  phenomenological  (or  kinetic)  coctfi  ieuts  which, in  the  subject  case  of 
absence  of  imposed  magnetic  field,  are  scil.tr  quantities.  The  Onsager's  re¬ 
ciprocity  relation  I,j  =  L^holds.  The  gradients  of  the  intensive  variables 
appearing  in  equations  (  l)  are  the  generalized  lories,  or  affinities. 

Different  alternative  forms,  each  suitable  lor  a  number  of  specific  uses, 
can  be  given  to  the  system  (4)  by  performing  suitable  linear  transformations 
of  the  affinities  and/or  the  fluxes,  according  to  a  well  defined  set  of  rules. 

If  Aj  indicates  the  direct  affinity  for  the  flux  Jj  and  if  a  prime  indicates  the 
transformed  quantities,  then  the  following  relations  hold:  (Ret.  li) 


A!  -  /•  P:uA 


ik 


p;Ljk 


L1  =  I,1  .  -  >  (3*1  L..  (3* '  • 

mi  nn  ■<-.  ■  ik  ik  mi 

J-  k 

where  P-^  is  the  reciprocal  of  the  element  Pik  in  the  square  matrix 
which  defines  the  linear  transformation  of  the  affinities. 
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We  shall  make  use  of  few  such  transformations  during  the  development 
of  the  analysis.  A  particular  transformation  is  however  very  instructive 
for  an  introductive  survey  of-the  nature  and  types  of  results  which  we 
will  get. 

We  notice  first  of  all.  that  the  local  thermodynamic  state;  ot  the 
mixture,  in  the  subject  case  of  thermal  equilibrium  among  the  con¬ 
stituents,  is  characterized,  as  follows  from  eq.  (  i)  when  accounting 
for  the  presence  of  the  electric  potential  6  ,  by  five  independent  variables. 
Go in>;  from  one  set  of  variables  to  another  is  accomplished  through  the 
"state  equations"  of  the  mixture.  It  proves  convenient  to  assume  the 
following  set  of  parameters  <£.  T,  p  (i=  1.  *i).  where  the  c.'s  are 
the  mass  concentrations  of  the  charged  components,  as  the  basic  set  of 
independent  thermodynamic  variables. 

The  gradients  of  these  variables  are  not  all  independent  since  there 
are  a  number  of  conservation  equations  to  he  satisfied  and  there  is,  in 
addition,  the  requirement  of  macroscopic  neutrality.  The  latter  amounts 
to  the  relation 

e,  V  i  ,  4  e,  V  i  ,  0  ( t>) 

between  the  two  concentration  gradients,  while  the  momentum  conservation 
requires,  under  the  subject  assumptions  of  macroscopic  neutrality  and 
absence  of  mass  motion,  that  the  gradient  of  p  vanish  identically  throughout 
the  field. 

One  is  thus  left  with  only  three  independent  gradients:  say  VT,  V  <j> 


and  V  cj 
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It  is  most  profitable  and  instructive  to  assume  as  new  generalized  forces 
these  three  independent  gradients  and  express,  accordingly,  the  dynamic 
equations  (eqs.  4)  in  terms  of  these  forces.  This  can  be  done  by  expressing, 
through  the  pertinent  state  equations,  the  gradients  appearing  in  eqs.  (4)  as 
linear  combinations  of  VT,V<j>  ,Vr(  and  then  applying  eqs.  (  5)  - 
As  a  result,  the  dynamical  equations  will  read: 


Jo  =  T 

t 

Lo, 

+ 

v<t> 

J,'  =  L'.o^T 

+ 

L.U 

Vc, 

+ 

l',2 

v<t) 

(7) 

j!  =  L,'0  v  r 

+ 

L>i 

Vc, 

4 

1.1  i 

V4> 

A'0  =VT 


A  2  = 


and  the  new  fluxes 


J.  will 


be  given  by: 


J 


i 

o 


(  S|  Jj  +  s,  J, 


+  s, 


Jl 


1 

<'(pi  -p  i ) 

T 

.  dci 

0  1*2 


J. 


i- 


X 


-P3 


I 


e,  0(g,  -pj) 

9(2 


(9) 


J 


i 

i 


j. 

T 


e,  J,  + 


T 
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where  s  are  the  specific  entropies  of  the  constituents,  and  Js  is  the  total 
entropy  flux: 

js  =  i!L  .  hJi  _  jLl£l  .  iLiii.  (io) 

s  j  T  T  T  '  ' 

The  quantity  j'Q  can  be  thought  of  as  the  "conduction"  flux  of  entropy. 

The  flux  Tj]  has  the  clear  physical  meaning  of  the  electric  current.  No 
direct  physical  interpretation  can  be  given  to  the;  flux  Jj  which  is  a  partic¬ 
ular  linear  combination  of  the  mass  fluxes  Jj  and  J.  .  The  thermodynamic 
derivate  appearing  in  eq.  (  9)  are  computed  holding  T  and  p  constant. 

The  l| j  '*>  are  the  new  phenomenological  coefficients  in  the  present 
context;  there  is  no  need  for  expressing  L?.  in  terms  of  t lie  original 
ones  L-j  . 

The  system  of  eqs.  (9)  ,  containing  the  three  truly  independent  gradients 
V  4>  ,  VT,  V  Cj  ,  is  in  the  appropriate  form  for  the  proper  application  of  the 
theorems  of  the  thermodynamics  of  irreversible  processes. 

Suppose  we  maintain,  from  outside  ,  a  constant  temperature  gradient 
on  the  boundaries  of  the  system. 

Eventually  a  steady  state  will  be  reached,  referred  to  as  a  stationary 
states  of  first  order,  in  which  the  direct  flux  is  the  only  flux  different 

from  zero.  System  (  9)  will  then  contain  only  four  "unknowns"  (the 
aforementioned  gradients  plus  )  so  that  three  of  them  can  be  expressed 

as  function  of  the  gradient  of  temperature.  One  thus  obtains  relations  of 
the  type: 

V<(>  =  u(  s)  V  T 

V  Cj  =  (  s)  V  T 


(ID 
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where  the  symbol  (  s  )  stands  for  "function  of  the  local  thermodynamic  state 
of  the  fluid". 

The  first  relation  gives  the  electric  field  induced  by  the  non -uniformity 
(under  the  restraint  of  constant  gradient  on  the  boundary) .  The  phenomeno 
logical  coefficient  u(s)  is  thus  one  of  the  required  electromagnetic  prop¬ 
erties"  of  the  medium.  The  second  relation  gives  the  gradient  of  concentra¬ 
tion  induced  by  the  temperature  gradient  in  the  presence  of  the  induced 
electric  field.  The  last  equation  relates  the  temperature  gradient  to  the 
energy  flux  into  the  system  (notice  that  =  .1,'  *  0  implies  absence  of  any 
mass  flux  throughout  the  system).  When  computed  at  the  boundaries,  it  can 
be  interpreted,  in  tile  light  of  an  overall  energy  balance,  as  the  energy  to  he 
supplied  from  the  "ambient"  to  mantain  the  induced  electric  field  (or,  what 
amounts  to  the  same,  the  induced  concentration  gradients)  . 

Equations  (9)  ,  as  they  stand,  can  be  used  for  order  of  magntitude 
analyses  of  the  effects  of  the  non-uniformity  herein  considered,  their  degree 
ot  accuracy  depending,  apart  from  the  approximations  involved  in  exploiting 
the  phenomenological  coefficients,  on  the  overall  dimensions  of  the  system. 
Obviously,  the  thermodynamic  state  of  the  medium  will  be  a  function  of  the 
point  so  that,  for  instance,  the  first  of  eqs.  (II)  should  be  written  as: 

V  <t»  (r)  =  <i  [*(r)J  VT(r)  [li) 

where  r  is  a  space  coordinate.  Thus  the  "complete"  solution  of  the  problem 
calls  for  the  determination  of  the  state  parameters  throughout  the  medium 
(in  order  to  answer,  for  instance  ,  the  question  of  what  is  the  value  of  V  T 
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throughout  the  system  when  a  given  non-uni*orm  temperature  distribution  is 
maintained  on  its  boundaries).  This  is,  obviously,  a  much  more  complex 
matter  which  must  be  handled  by  solving  the  energy  and  mass  equations 
JG=  ^i  =  ^  subject  to  suitable  conditions  on  the  boundaries  of  the 
medium. 

However,  order  of  magnitude  analyses  can  always  be  carried  out  by 
taking  suitable  "average"  values  of  the  state  parameters. 

Since,  in  this  case  J,  «  jJ  =0  implies  also  J,  =  J2  =0  .this  type  of 
information  can  also  be  obtained  directly  from  the  system  (4)  by  setting  in 
it  J,  =  J,  =  0.  This  simpler  procedure  will  indeed  be  followed  in  the  next  part. 

Consider  now  the  other  case  in  which  a  given  non-uniform  distribution  of 
electric  potential  O  is  maintained  on  the  boundaries.  As  before,  the  principle 
of  stationary  states  implies  that  the  only  flux  different  from  zero  is  the 
direct  flux  J,'  that  is  the  electrical  current  I.  Notice  that  neither  the 
energy  flux  nor  the  mass  flux  is  singularly  equal  to  zero:  what  vanishes  is 
rather  two  linear  combinations  of  them,  as  expressed  in  terms  of  j'Q  and 
J2  .  This  clearly  shows  why  the  principle  of  stationary  states  could  have  not 
been  applied  to  the  dynamical  equations  in  the  form  given  by  equations  (  4)  . 
System  (  9)  contains  once  again  only  four  unknowns  so  that  it  can  be  solved  to 
obtain  the  following  relations: 

I  -  CT  (  s)  V* 

VT  -  (  s)  V<t> 

Vc,  =  (s)  V<t> 


(13) 


■ 


40 

The  first  equation  is  of  foremost  interest  for  the  present  analysis.  It  yields 
the  steady  state  "electrical  conductivity "  of  the  non-uniform  medium  subject 
to  the  conditions  of  a  given  distribution  of  electrical  potential  on  its  boundaries. 
In  the  prec.eeding  case  an  energy  exchange  between  the  medium  and  its  sur¬ 
roundings  was  necessary  to  maintain,  in  the  stationary  case,  the  induced 
electric  field.  Now  an  exchange  of  electric  charges  between  the  medium  and 
its  surroundings  is  necessary,  its  amount  being  determined  by  eq.  (  13) 
evaluated  along  the  boundaries  of  the  region.  Once  again,  complete  quanti¬ 
tative  solution  of  the  problem  calls  for  the  solution  of  the  mass  and  energy 
conservation  equations  to  determine  the  actual  distribution  of  the  state  para¬ 
meters  throughout  the  system  considered. 

This  is  the  type  of  information  which  can  be  obtained  from  the  present 
analysis.  The  remainder  of  the  body  of  the  chapter  is  devoted  to  obtaining 
explicit  expressions  for  the  coefficients  u  (  s)  and  IT  ( s) . 


I. 

I 

[ 

[ 


L  -  Deter  ruination  of  the  coefficient  n  , 

From  the  last  section  it  follows  that  to  determine  u  one  must  consider 
the  stationary  state  for  which  a  given  nonuniform  temperature  distribution 
is  kept  constant  on  the  boundaries  of  a  system  and  the  mass  fluxes  J,  and 
j2  vanish  throughout. 

It  here  proves  convenient  to  adopt  the  following  form  of  the  dynamical 
equations: 


\r-,  Qi  Ji 


i=  l 


+  (Leo  -  L  jo  Qj)  Au 


(14) 


ji  =  Lik  <  V  A>  +  QkA! 

k=  l 


I 
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where  the  affinities  A.  are  now  defined  as: 

1 


A 


u 


T 


V  T 


A  i  =  -  c,  V  -  TV  (J^_) 


A , 


e,  V  4-  -  T  V 


up 


A,  =  -  TV  (H_l) 


(15) 


and  where  the  quantities  Q.(  i=  1.2)  are  the  heats  of  transfer  (Ref.  19) 
defined  in  terms  of  the  coefficients  L..  appearing  in  eq.  (4),  as: 


l^io  l^o  ”  l'ii  ^io 

G,  - - 2 - 

1  ..  i  -  L,- 


(16) 


‘Ml  1 J- o  "  I'lO  1-i  1 
^-1 1  1-ii  "  1-u 


Express ing  the  gradients  of  p.  in  terms  of  the  independent  gradients 
V  r,  v4>,Vc,  .and  proceeding  as  indicated  in  tin-  1  »st  section  it  can  be  shown 
that  the  coefficient  a  has  the  following  expression: 


(li  (  b i  /ei  “  ai  / ei  )  +  <i>  (  ij/ 1 ,  -  b, /cj  ) 


(17) 
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with: 

ai  =  9(  ^  *  p ,  )  /  9c, 

-  H-i)  /  9c, 

(18) 

®i*  n 

q:  -  —  (U;  -  >>i  ‘  M 
T 

whoro  Vj  and  h{  are  the  specific  volume  and  enthalpy  of  the  i-th 
constituent ,  respectively. 

It  appears  that  the  aj.  1)..  m;  are  completely  determined  once  the 

pertinent  "state  equations'  (Ref.  21)  of  the  constituents  are  given.  The 

quantities  cjj  depend,  in  addition,  upon  the  phenomenological  coefficients 

L  through  the  heats  of  transfer  Qj  . 

U 

Thus,  to  make  any  use  of  cq.  (17)  one  must  make  some  assumptions  as 
to  the  state  equations  of  the  plasma  constituents  and  must  determine,  either 
experimentally  or  by  means  of  statistical  mechanics,  the  phenomenological 

coeffic  ients  L,  j  . 

The  expressions  for  a*  .  b.  and  m.  are  herein  derived  on  the  assumption 
that  the  plasma  constituents  are  perfect  gases.  1' he  quantities  q,  will  be 
dealt  with  in  the  next  paragraph. 
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If  one  makes  the  perfect  gas  assumption,  the  following  relations  hold 
(  ref.  23)  . 


Ri  = 


RT 


Ill; 


4>:  (T)  +  l 


"  fa  ^f) 


(19) 


wherein  (*>.  (T)  is  an  arbitrary  function  of  the  temperature,  irij  is  the 
molar  mass  of  the  it h  component,  m  is  the  molar  mass  of  the  ith  com¬ 
ponent,  m  is  the  molar  mass  of  the  mixture,  h0  a  reference  specific 
enthalpy,  and  R  the  molar  gas  constant. 

With  these  expressions  one  has: 
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mi 


In  the  case  of  an  imperfect  Lorentz 's  gas  (  i.  e.  /mi«l  ;  N,/N  «1; 

1  3 

)  these  formulas  are  greatly  simplified  and  one  obtains: 


a,  RTN  m 

e,  t ,  m, 


b,  =  a,  =  U 


b,  RTN  m 

t,  *2  m> 


nii 


RT 

pm 


(Zl) 


where  <j  =  — Vr  are  the  electric  charges  per  unit  volume. 
mi*S 

The  explicit  expression  for  the  coefticients  u  .  is  lound  by  substitutiing 
eqs.  (  ZO)  into  eq.  (17). 

In  the  case  of  imperfect  Lorentz's  gas  one  obtains: 


a 


q. 


Ntn  JjHIL. 


«i'"i 


( 2  In2i ,  t,tn,c2 


(ZZ) 


The  formulae  so  lar  developed  are  all  that  one  can  do  without  specifying 
the  nature  of  the  phenomenological  coefficients  Ljj.  Any  order  of  magnitude 
analysis  can  be  furthered  only  after  having  obtained  suitable  expressions  for 
these  coefficients  or,  what  amounts  to  the  same,  for  the  quantities  qp 
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This  is  considered  in  the  next  part,  where  it  will  be  shown  how  the 
L-'s  can  all  be  expressed  in  terms  of  transport  coefficients. 

5.  Expressions  of  the  pertinent  phenomenological  coefticients  in 
terms  of  transport  coefficients 

To  express  the  coeff ic ient s  Ljj  in  terms  oi  transport  coefficients  the 
following  consideration  is  essential. 

In  the  subject  case  we  have  two  independent  mass  fluxes  and  one  energy 
flux.  The  number  of  independent  phenomenological  coefficients  is  therefore 
six,  accounting  for  the  three  Onsager  reciprocity  relations.  It  follows  that 
the  behavior  of  the  system  is  com pletely  characterized  by  only  six  indepen¬ 
dent  coefficients.  One  such  a  set  is  that  of  the  three  binary  diffusion  co¬ 
efficients  D,z,  D  ,,  D,j  (  i.e.  ,  coefficients  of  diffusion  of  electron  -  ion , 
electro-molecule  and  ion-molecule  mixtures)  plus  the  two  thermal  diffusion 
coefficients  D1  ,  1)2  (i.e.  the  coefficients  of  thermal  transport  due  to  dif¬ 
fusion  of  the  electron  gas  and  the  ion  gas  into  the  neutral-ion  and  neutral- 
electron  gas,  respectively)  and  the  heat  conduction  coetlicient. 

It  is  then  natural  to  think  that  a  suitable  linear  transformation  of  fluxes 
and  forces  will  make  it  possible  to  express  everything  in  terms  of  the  above 


transport  coefficients. 


45 


This  is  considered  in  the  next  part,  where  it  will  be  shown  how  the 
L.  .'s  can  all  be  expressed  in  terms  of  transport  coefficients. 


3.  Expressions  of  the  pertinent  phenomenological  coefficients  in 
terms  of  transport  coefficients 

To  express  the  coeffic ients  Ljj  in  terms  of  transport  coefficients  the 
following  consideration  is  essential. 

In  the  subject  case  we  have  two  independent  mass  fluxes  and  one  energy 
flux.  The  number  of  inde pendent  phenomenological  coefficients  is  therefoi  e 
six,  accounting  for  the  three  Onsager  reciprocity  relations.  It  follows  that 
the  behavior  of  the  system  is  completely  characterized  by  only  six  indepen¬ 
dent  coefficients.  One  such  a  set  is  that  of  the  three  binary  diffusion  co¬ 
efficients  D  2.  D  , ,  n^j  (  i.e.  .  coefficients  of  diffusion  of  electron-ion, 
electro-molecule  and  ion-molecule  mixtures)  plus  the  two  thermal  diffusion 
coefficients  D,1  ,  1)2  (i.e.  the  coefficients  of  thermal  transport  due  to  dif¬ 
fusion  of  the  electron  gas  and  the  ion  gas  into  the  neutral-ion  and  neutral- 
electron  gas,  respectively)  and  the  heat  conduction  coefficient. 

It  is  then  natural  to  think  that  a  suitable  linear  transformation  of  fluxes 
and  forces  will  make  it  possible  to  express  everything  in  terms  of  the  above 


transport  coefficients. 


4b 


This  is  indeed  so.  In  fact,  by  defining  new  generalized  forces  A-  as: 
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4-  h; 


A, 


ii 
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(i  =  l,  2,  2  ) 


(22) 


we  find,  according  to  eqs.  (.24),  the  following  expressions  for  the  new 
fluxes: 


(i  =  1.  2,  2  ) 


and  the  new  phenomenological  coefficients  L!  : 
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These  new  coefficients  L"  are  just  those  used  in  reference  (  1)  and  it  is 

U 

therein  shown  how  they  are  related  to  the  binary  diffusion  coefficients  Djj. 
to  the  thermal  diffusion  coefficients  D T  (  i=  1  ,  2)  and  to  the  thermal  con¬ 
ductivity  coefficient  X.  . 

The  ultimate  goal  of  expressing  the  qj's  in  terms  of  transport  coefficients 
is  thus  achieved  by  performing  the  necessary  substitutions.  One  gets 
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and  by  using  the  expressions  for  L;,  reported  in  Ref.  (  1 ) : 
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and: 


q»  =  * 
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4.  -  Determination  of  the  electric  conductivity  of  a  plasma 
As  shown  in  the  proceeding  parts  the  determination  of  the  coefficient  a 
in  terms  of  the  transport  coefficients  of  the  mixture  involves  essentially  two 
steps.  First  the  dynamical  equations  are  written  in  the  form  which  is  appro¬ 
priate  for  the  subsequent  application  of  the  principle  of  stationary  state  which 

i 

yields  the  expression  of  ci  in  terms  of  some  set  of  kinetic  coefficients  Ljj 
and  some  thermodynamics  derivatives.  To  express  the  latter  ones  some 
assumptions  must  be  made  as  to  the  nature  of  the  state  equations  of  the  con¬ 
stituent  gasses  (in  particular,  the  perfect  gas  hypothesis  has  been  made). 

In  the  second  step,  suitable  linear  transformations  are  introduced  which 

*  I 

related  the  kinetic  coefficients  L-(  to  a  new  set  of  coefficients  Ljj  which,  as 
shown  in  Ref.  (  1)  ,  is  directly  expressible  in  terms  of  the  transport  coefficients. 

This  same  procedure  must  be  followed  to  determine  the  electrical  conduc¬ 
tivity  0  ,  and  can  be  reported  now  in  a  condensed  form,  by  skipping  all 

details. 

The  appropriate  form  of  the  dynamic  equations  is  that  given  in  eqs.  (7)  . 

By  definition,  the  electric  conductivity  given  by  the  ratio  i/v*  between  the 
electric  current  density  and  the  electric  tield. 

This  ratio  depends  on  the  particular  assumptions  about  the  order  of  the 
stationary  state  in  which  it  is  evaluated  and  thus  to  a  certain  extent  on  the 
number  of  "constraints"  which  arc  imposed  at  the  boundary  of  the  system. 

It  acquires  a  definite  meaning,  that  is,  it  can  be  given  a  unique  expression 
in  terms  of  the  state  parameters,  only  in  a  particular  first  order  stationary 


state. 
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Indeed,  in  the  subject  case  of  a  ternary  plasma,  one  could  define  a 
third  order  stationary  state  in  which  a  given  non-uniform  electric  potential 
distribution  and  a  uniform  temperature  and  electron  concentrations  are 
maintained  on  the  boundaries  of  the  system.  In  such  a  state,  however,  the 
plasma  would  certainly  not  be  uniform  throughout  the  region  (i.e.VTand  Vc, 
would  not  be  identically  zero)  and  the  ratio  between  I  and  V  4'  would  also  be 
a  function  of  the  other  gradients  (and/or  fluxes)  present  in  the  system. 

Similar  considerations  hold  for  a  second  order  stationary  state  in  which  a  given 
non-uniform  electric  potential  distribution  and  a  uniform  temperature  (or 
electron  concentration)  are  maintained  on  the  boundary  of  the  system. 

Once  again,  the  gradient  of  T  (or  of  c,)  would  not  vanish  identically 
throughout  the  system,  two  independent  fluxes  would  be  different  from  zero 
and  no  unique  definition  for  the  electric  conductivity  could  be  given  (i.e.  the 
ratio  1/  V  <*  would  not  depend  only  on  the  state  of  the  system  but  also  on  some 
gradients  of  the  state  parameters). 

Thus  the  electric  conductivity  must  be  defined  in  a  first  order  stationary 
state  tor  which  a  given  non  uniform  distribution  of  electric  potent ive  is 
maintained  on  the  boundaries  of  the  system.  In  this  case  the  electric  current 
will  be  the  only  non -vanishing  flux  and  equations  (7)  can  be  solved  for  the 

ratio  0  :  (l/V6)  giving: 
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It  is  of  interest  to  notice  that  in  a  fully  ionized  neutral  plasma  (i.e.  in  a 
system  with  one  less  degree  of  freedom)  in  which  the  electron  concentration 
gradient  V  is  no  longer  an  independent  gradient,  the  expression  for  0  in 
the  first  order  stationary  state  becomes  by  formally  putting  L(V  -  0  in  Eq.(7): 


d  =  L 


L«2 


(30) 


'oo 


This  relation,  when  expressed  in  terms  of  the  only  binary  diffusion  coeffic¬ 
ient  existing  in  this  case,  reduces,  as  it  is  easy  to  verify,  to  the  well 
known  Einstein-type  relation. 

To  obtain  U"  in  terms  of  the  set  of  transport  coefficients,  one  only 
needs  a  set  of  relations  between  the  Ljj  and  the  's,  which,  as  said 
before,  are  related  to  the  transport  coefficients  as  shown  in  Ref.  (  1)  The 
subject  relationships  are: 
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where  all  the  symbols  have  the  already  defined  meaning. 
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Z.  1 .  -  The  transport  coefficients  for  a  plasma. 

In  the  proceeding  section  the  general  expressions  for  the  coefficients 

T  T 

a  and  O’  interims  of  the  transport  coefficients  Dj>,  Dj  ^ ,  D^,  ,  D^.X. 

Iiave  been  derived. 

The  next  step  toward  the  ultimate  goal  of  expressing  a  and  0  ,  to 
within  a  given  degree  of  approximation,  in  terms  of  the  thermodynamic 
parameters  of  the  mixture  and  of  the  physical  charateristics  of  the  con¬ 
stituent  gases  requires  then  the  evaluation  of  the  above  mentioned  transport 
coefficients  in  terms  of  these  parameters. 

This  will  be  done  in  this  section. 

We  recall  some  basic  facts  on  the  determination  of  the  transport  co¬ 
efficients. 

For  conditions  not  too  far  from  equilibrium  the  determination  of  the 
transport  coefficients  hinges  on  the*  evaluation  of  the  following  integrals 
(see  Refs.  1  or  6  details): 

I)  [’lie  deflection  angle: 
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and  where:  mjj  is  the  reduced  mass,  defined  as  m^-  whe rc  mj 

and  mj  are  the  masses  of  the  particules  of  species  i  and  j;  g  is  the  initial 
relative  velocity  of  the  two  particles;  b  is  the  impact  parameter;  r  the 
distance  between  the  mass  centers  of  the  two  colliding  particles;  <?  (  r)  is  the 
interaction  potential;  rm  the  distance  ot  closest  approach. 

Eq.  (3  2)  is  valid  for  elastic  collisions  between  particles  assumed  to  be 

centers  of  force  fields  with  spherical  symmetry. 

The  potentials  which  we  shall  deal  with  are  all  expressible  as: 

Mr)  :  (J) 


where  «  and  0  depend  on  the  particular  'molecular  model  assumed.  In 

particular,  0  corresponds  to  the  molecular  diameter  in  the  case  ot  the 
rigid  sphere  model. 
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•  re  the  subscripts  r.  s.  indicate  values  for  the  rigid  sphere  model. 
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with  0  defined  as  before. 
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3)  The  non  dimensional  collision  integrals; 
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Notice  that,  by  definiton,  all  starred  quantities  are  equal  to  one  for  the 
rigid  sphere  model. 

Any  transport  coetlicienl  is  expressible  in  terms  of  the  above  quantities. 
Indeed,  indicating  by  the  subscripts  i,j,  quantities  related  to  collisions 
between  particles  of  the  ith  and  jth  spei  ies  one  has; 

«•»  »ij  -  (  diffusion  coefficient  for  the  ith  and  jth  species  in  a  multi- 
component  gas  mixture)  is  related  to  the  binaty  diffusion  coeffici 
ents  D- j  of  each  pair  of  species  (see  eq.  8-Z.49,  in  Ref.  (1)). 
These  relations  involve  in  addition,  the  concentrations  of  the 
single  species,  and  the  molecular  masses. 

The  coefficient  D.j,  in  turn,  is  expressible  as; 

D  -  ~  J_  A  mi  i  KT 
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is  the  collision  diameter. 


55 


b)  The  coefficients  of  thermal  diffusion,  D-  ,  can  be  evaluated  in  terms 
of  the  coefficients  of  binary  molecular  diffusion  and  of  the  coefficients 
of  thermal  conductivity  for  each  species  ( K^)  and  for  each  pair  of  species 
(  )  .  The  pertinent  explicit  expressions  are  too  cumbersome  to  be 

reported  here;  they  can  be  found  in  Ref.(  1)  page  543. 

It  is  important,  however,  to  notice  that  they  involve  the  integrals 


and  S>*,  (  s  =  t,  3,  3). 


ij  >J 

c)  The  coefficients  of  thermal  conductiv  ity  for  the  ith  species  can  be 


evaluated  as: 
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This  expression  is  valid  for  monatomic  gases;  Euckeri's  corrections  are  to  be 
added  for  polyatomic  gases. 

The  coefficient  of  thermal  conductivity  for  a  binary  mixture  is  given  by: 


P?Jj. 


(35) 


where  A;j  is  a  combination  of  the  collision  integrals  .  Its  val  le  is  very 

close  to  one. 

Table  I  summarizes  the  fundamental  dependence  ol  the  discussed  trans¬ 
port  coefficients  for  multicomponent  gaseous  mixtures  (first  column)  on  the 
related  binary  transport  coefficients  (column  two)  and  on  the  related  collision 
integrals  (column  three). 

The  next  step  is  a  discussion  of  the  actual  evaluation  of  the  transport 
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coefficients  for  the  particular  potentials  related  to  the  interactions  between 
the  particles  present  in  a  plasma, 
a)  Interaction  electron-ions 

The  force  between  two  charged  particles  is  coulombian. 

The  interaction  potential  is  given  by  the  equation: 


Ze,2 

*(«■)  =  -  ~ 


(36) 


where  Z  indicates  the  number  of  electric  charges  carried  by  the  ion  .and 
the  sign  takes  into  account  the  attractive  nature  of  the  force  for  dis¬ 
tances  greater  than  the  threshold  ot  actions  oi  quantic  type,  for  which  the 
force  is  repulsive. 

With  the  potential  given  by  eq,  (  3(/i  one  faces  the  well  known  fact  that 
the  collision  cross  section,  as  defined  by  eq.  (  32)  assumes  an  infinite 
value. 

Tills  problem  lias  been  studied  extensively,  and  the  difficulty  eliminated 
by  a  'cut  off"  of  the  upper  limit  of  the  integral.  The  cut  off  distance  must  be 
the  smaller  of  the  mean  distance  between  tin  gas  particles  D=  n  and  the 

Debye  shielding  length: 

>  1/2 

-  (KT/4  *  n,  e,  ) 

Consider  at  first,  as  pertinent,  the  collisions  for  which  the  impact 
distance  b  is  less  than  D. 

One  evaluates  eq.  (  32)  wherein  the  upper  limit  becomes  b=  D. 

Substituting  for  the  expression  ^  m  its  mean  value  2kT  and 

neglecting  the  small  variations  of  the  kinetic  energy  of  the  particles  within 
the  sphere  of  radius  D,  one  obtains: 
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<*"  <«>  *  7  (4^  '•> 


1  4  4 


2KTD  2 


V  Zc,Vj 


where,  for  mixture  of  electrons  -  ions  it  is: 

-*/, 


I)  =  n,  +  nz 

J 


(37) 


(38) 


The  corresponding  binary  diffusion  coefficient  is  obtained  by  sub¬ 
stituting  eq.  (  37)  into  eq.  (  33)  which  results  in: 


l< 

r,  . 

3  ( lrh:  KT 

;  i 1 KT  f  j\ 

<-  <n 

2KTD\* 

1  +4  j 

\ 

«  V  n  ) 

\  ‘  i  j  p , 

1  * 

'  e,  ez  J 

, 

L 

_  a 

J 

In  the  second  case,  when  the  upper  limit  is  put  equal  to  1D  tlie  col¬ 
lision  cross  section  will  be: 


(39) 


Q 


(i) 


2  ire. 


('/.  ki  r 


In  A 


whe  re: 


A  = 


L  1 


»  fjdrl'i 

-  e,  ’  tj  n, 


(10) 


1'he  corresponding  expression  for  Dj  ^  is  readily  obtained  by  substituting 
eq.  (  40)  into  eq.  (3  3)  and  is: 


D 


ii 


(41) 


More  recently  a  new  method  to  evaluate  the  coefficient  Dj  ^  has  been 
proposed  (Ref.  12)  which  does  not  use  the  cut  off  for  evaluating  the  collision 
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cross  section.  Such  a  method  furnishes,  with  a  more  complicated  process, 
an  approximation  higher  than  that  necessary  for  our  purposes,  and  it  is 
therefore  not  used  here. 

k)  Interaction  of  electrons -neutral  molecules  and  ions-neutral 
molecules 

The  interaction  between  charged  particles  (  electrons  or  ions)  and  the 
neutral  molecules  can  be  treated  by  the  following  procedure  (Ref.  i) 

The  charged  particle  produces  an  electric  field  that  induces  an  electric- 
dipole  into  tiie  neutral  molecule,  which  undergoes  a  separation  of  electric 
charges.  The  induced  electric  dipole  produces  an  external  potential  given 
by: 

V  =  +  a  —  cos  <R>  (i  =  l,  Z) 

r4 

where  a  is  the  polarizability  of  the  molecule,  Oj  the  inducing  charge,  and 
©  the  position  angle  referred  to  the  axis  of  the  dipole. 

The  interaction  potential  will  then  be  written  as: 


<?  (  r) 


HZ) 


The  negative  sign  indicates  tiiat  the  forces  are  attractive  up  to  the  distance 
,  quantic  threshold.  It  should  be  noted  that  the  energy  of  interaction 
expressed  by  ( -1Z)  is  only  a  first  approximation  to  more  complicated  express¬ 
ions,  because  it  neglects,  for  instance,  tiie  London's  dispersion  energy, 
related  to  the  mutual  action  between  the  molecular  dipole  and  the  dipole 
induced  by  this  into  the  ion.  Such  a  dispersion  energy,  proportional  to  the 


\ 
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inverse  sixth. power  of  the  distance,  can  amount  to  25%  of  the  interaction 
energy  given  by  (42)  for  distances  between  the  particles  of  the  order  of  one 
collision  diameter  (  Ref.  1 .  ch.  13). 

The  polarizability  of  the  atoms  and  of  the  molecules  has  been  studied 
theoretically  and  experimentally;  its  values  vary,  tor  polyatomic  molecules, 
with  the  direction  of  the  inducing  electric  field.  It  is  always  possible  however 
to  introduce  a  mean  value  taken  over  all  directions.  This  value  increases  with 
the  molecular  complexity;  for  instance,  the  mean  polarizability  varies  from 
<•  ’  x  10  cm  to  103,2  x  10“  cm3  for  the  molecule  of  hydrogen  and  the 
molecule  of  henzone,  respectively.  Tables,  given  in  Ref.  (  1)  furnish  the  more 
important  data. 

l  or  a  potential  expressed  by  eq.  (42)  which  is  nothing  but  a  particular 
form  of  the  Sutherland's  model  with  y=  4.  one  notices,  from  the  tab.  4  of 
Ret.(  1)  that  tor  not  too  high  values  ot  the  particle  energy,  the  collision 
cross  sections  become  very  close  to  the  collision  cross  sections  tor  rigid 
spheres.  For  first  approximation  calculations  and  for  orders  of  magnitude 
evaluations  and  in  view  of  the  uncertainty  as  to  the  value  of  „  ,  we  ran 
the  retore  use  the  results  obtained  for  the  simple  potential  of  rigid  spheres 
without  too  large  errors. 

We  write  therefore; 
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1 1 
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i  i 


-  3 

1 

KT  1 

8 

V 

1  t  m  ( (  n 

(4  3) 

=  3 

l 

[  KT*  l 

8 

<  \i 

11 

(44) 

t.0 


where: 


“u  ’  4-  <0 1  ’  °J » 

is  the  collision  diameter. 

An  alternate  procedure  for  the  determination  of  D}-}(  i=  1,2), 
based  on  experimental  data,  can  be  used. 

Recently,  experimental  evaluations  of  the  total  elastic  collision  cross 
sections  were  made  for  the  ions-neutral  molecules  and  the  electrons -neutral 
molecules  mixtures. 

The  total  elastic  collision  cross  sections  are  defined  by: 


(°) 

Q  =  2  u  l  b  ti  b 


(45) 


and  one  can  write: 


Q(°  =  Q' 


(46) 


where  (1-  cos  )  indicates  the  mean  val  ic  of  the  term  (  I -cos  f  )  taken 
over  all  angles  oi  collision.  Such  a  mean  val  le,  lor  collision  angles  equally 
probable,  is  equal  to  one.  This  occurs  for  not  too  elevated  values  of  the  energy 
of  the  particles. 

With  this  limitation  and  with  not  too  large  errors,  one  can  then  assume 

/  |  ,  (  o ) 

for  Q  the  values  experimentally  obtained  tor  Q  and  hence  compute 

□  i  3  and  by  substituting  these  values  in  eqs.  (  43-44). 

A  more  extended  discussion  of  this  procedure  can  be  found  in  Ref.  9. 
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c)  Interaction  between  neutral  molecules 
The  interaction  potential  is  expressed  by: 

<t>  (  r)  =  -  (47) 

r 

since  the  two  molecules  behave  as  two  mutually  inducing  dipoles.  This 

potential,  corresponding  to  the  Sutherland's  model,  was  studied  by 

(es>* 

Kotani  (Ret.  4)  and  the  values  of  .i  are  available*  in  tabulated  form. 

The  evaluation  of  the  transport  coefficients  for  this  model  is  thus 
reduced  to  the  substitution  ot  these  tabulated  values  into  the  previously 
described  formulae. 

The  results  thus  far  described  arc*  summarized  in  Table  II.  For  each 
of  the  three  types  of  interactions  the  cor responding  expressions  for  the 
potential  are  reported  together  with  the  available  references  for  the  theoretical 
and/or  experimental  evaluation  of  the  collision  cross  sections. 

i..  -  Analysis  of  the  orders  of  magnitude*  of  the  transport  coefficients 
The*  determination  of  the  transport  coeff ic ients .  as  described  above, 
will  now  lie*  specialized  to  the  case  of  a  weakly  ionized,  imperfect  Lorentz's 
gas  for  which,  ns  said,  the  following  inequalities  hold  by  definition: 


nil  ~  mi 


«  1; 


--«1 


11(2  in3  ‘*3  *>j 

It  will  be  found  that  in  such  case*  a  number  of  terms  and/or  contributions  can 

be  neglected  because  of  smaller  order.  To  lend  rigorousness  to  such  an  order 

of  magnitude  analysis,  we  shall  assume  ns  a  reference  "order  of  smallness" 

1/ 

the  square  root  of  the  mass  ratio  6  =  (  nij/iiq)  ,  For  the  neutral  particles 

.  2 

we  shall  be  interested  in,  6  is  of  the  order  of  1  0 


Let  us  begin  with  the  biliary  diffusion  coefficients.  It  is  pointed  out  that 
evaluation  of  their  comparative  orders  of  magnitude  will  be  made  for  given 
temperature  and  number  density  of  the  plasma. 

Eqs.  (  39-41-43-44)  indicate  that  for  given  n  and  T  the  order  of  magnitude 
of  the  coefficients  of  binary  molecular  diffusion  depend  on  both  the  ratio  o*  the 
colliding  molecular  masses  .and  the  ratio  of  the  related  collision  cross  sections. 

The  reduced  masses  in  the  present  approximation  are  given  by: 

nij,  m, 

m.  m, 

m2i 

The  orders  of  magnitude  of  the  ratios  between  the  collision  cross  sections  are 
immediately  eval  mted  for  any  energy  level  (this  being  the  only  parameter  on 
which  the  Qj  (  depend)  .  As  discussed  in  Section  II,  can  be  obtained  theor¬ 

etically,  while  either  rigid  sphere  model  formulae  or  experimental  results, 

|  reported,  tor  instance  ,  in  Hot's.  (  I  1 , 9 ,  3  ,  L  1  ,  I  8)  j  can  he  used  tor  Qj  j  • 

It  is  important  to  observe  that  the  ratios  between  the  collision  cross 
sections,  to  a  first  approximation,  do  not  depend  on  the  state  parameters. 
Indeed,  consider,  for  instance,  the  ratio  *n  a  plasma,  when 

n i  —  n^  =  106  part.  /cm\ 

In  the  range  1  ,  000<  T<  10,000°K,  such  a  ratio  varies  between  1.68x10“ 
and  1.10x10  .  For  n,=  n  =  10  part/cm  ,  the  ratio  is  roughly  doubled.  The 
above  reported  values  for  the  electron  number  density  are  limiting  values  for 
the  range  of  variation  of  n^,  corresponding  to  the  weakly  ionized  gases  of  the 
ionospheric  F  layer  and  of  the  electric  strong  discharges,  respectively. 
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As  a  consequence: 


Q2I  /  0\  2  = 

0  | 

K)-2]  s 

0  [  J 

Qa  !  Qi  j  = 

0 

['°2  ]  • 

ot"1] 

(48) 

u,,  /  Q,2s 

0 

These  particular  values  have  been  obtained  for  a  plasma  whose  con¬ 
stituents  are  N,,  N  }  ,  e  ,  and  at  1,000°K,  itj  n;  10^  part/cm\  The 
orders  of  magnitude  ot  these  ratios  do  not  essentially  vary  as  functions  of 
the  state  parameters  or.  for  the  usual  gases,  of  the  chemical  nature  of  the 
pa  rticlcs. 

Thus,  we  obtain,  for  the  diffusion  coefficients: 

D.  2  /  D2j  =  0  (  1 J 

D,  ,/  D„  =  0  [  6*^J  (4  9) 

I),,/  D,  3  =  U  jy] 

,  r  t 

(  b)  The  evaluation  ol  Dj  ,  and  ot  their  ratio  now  constitutes  the 

second  step. 

As  mentioned  previously,  the  coefficients  of  thermal  diffusion  for  each 
consitiuent  species  of  the  mixture  is  a  function  of  the  coefficients  of  binary 
molecular  diffusion  D-  and  of  the  coefficients  of  thermal  conductivity 

\ j  •  i  \ )  • 

Since: 

Qi  / Qt  =  o  y] 

Q;  / Qi  =  o  y] ' 

/  i  1  1 

Qi  /Q  •  5  0  L1  -j 


(50) 
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one  lias: 


\,/\ , 

=  0 

M 

VjAj 

=  0 

Xi  A> 

=  (J 

•  -*  T  T 

The  various  terms  in  the  equations  for  I)  j  and  can  lie  noticeably  simplified 
from  (  49)  and  (  51 )  . 

Neglecting  all  the  necessary  long  manipulations,  we  only  state  here  the 
pertinent  results  as  follows: 

1>,T  /  n?  =  D  (53 

liT  =  11  mn  7n  n,,  (5  i 

4  n,  v- 

who  re : 

ff[  =  M_.  ( I  i  M ,/  )  j (11  M,/M,4  3  5/4  Mj/Mj  +  13) 

For  convenience,  a  list  of  the  relative  orders  of  magnitude  of  all  the 


quantities  involved  is  reported  in  TABLF.  Ill, 


1 

I 

I 

1 

I 
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SECTION  3 


Oil  the  basis  of  the  results  of  the  last  sections,  one  can  examine  more 
accurately  the  phenomenological  coefficients  describing  the  electrical  prop¬ 
erties  of  a  non-uniform  plasma.  Namely,  at  this  point  it  is  possible  to 
evaluate  their  expressions  in  terms  of  the  ther motlynamic  state  parameters 
of  the  plasma. 


1 )  T he  Coefficient  n  ■ 

In  section  I  we  have  obtained  for  u  ,  lor  the  imperfect  Lorentz's  gas, 
the  expression: 


ru  hi  mi  +  R  i  A  i  m> 

n  =  - - 

tj  /<.,  m,  m  ,n,  +  <> /<  in>i  m2n2 


(54) 


Since  it  is  possible  to  evaluate  for  <|j,  and  q^,  a  relative  order  of 
magnitude  q^/qt-  0  by  analyzing  their  expression  (  lb-17),  one  obtains 

as  a  consequence: 


<12  f  ,  n,  n ,  m  2 


where  q^  assumes  now  the  following  form: 

11 1  D,,' 


^  a 


T 

A  K  D. 
m,  n,  Dzi 


l  + 


rb  ^i; 


Thus: 


n.  D»  ^ 

=  -  KP  Am  F  l  +  —  --ii 

e  V  “)  ^ 


(55) 


(5b) 


(57) 
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where: 


f  =  /  —  jnl  fa-U:  fa  Sl  +  fa 

\4  m,m2  n3  /  n2  cz  ei 


The  coefficient  a  is  dimensionally  expressed  in  volt /cm / °  K/cm  and 
represents  a  measure  of  the  electric  field  induced  by  a  unit  gradient  of 
temperature  in  a  plasma. 

It  has,  for  an  ionized  gas,  the  same  significance  of  the  Seebcck's  coefficient, 
usually  treated  for  solid  substances.  It  is  recalled  that  the  Seebeck  s 
coefficient  S  for  solid  conductors  is  defined  in  the  same  way  as  the  coefficient 
n  .  Indeed  S  is  the  ratio  between  the  gradients  of  electric  potential  and 
temperature  in  a  first  order  stationary  state  in  which  no  electric  current 
flows.  (See  Tef.23,  for  instance)j  . 

Examining  the  expression  for  a  one  observes  that  the  only  fundamental 
parameter  upon  which  it  substantially  depends  is  the  degree  of  ionization. 

Indeed,  for  a  weakly  ionized  gaseous  mixture,  one  can  suppose  that: 


m 

Tm 


;  are  only  slightly  different  from  unity;  and 


n2  ■  1  e2 

the  coefficient  K  AntF  is  a  constant  depending  on  the  mean  molecular  mass 
m  =  «  j  mjnj/n  only. 

The  ratio  between  the  coefficients  of  the  binary  molecular  diffusion 
D  /Di  as  we  have  shown,  does  not  depend  on  any  state  parameter.  Tin: 
only  parameter  that,  as  a  consequence,  affects  u  is  the  degree  ot  ionization. 
An  indicative  value  obtained  for  a  in  a  particular  interesting  case  follows. 
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2-  Value  of  a  obtained  for  the  ionosphere  F  layer 

Considering  a  plasma,  whose  characteristics  are  those  of  the  air  in  the 

ionspheric  F  layer,  we  have: 

T  =  l,  000°  K 


ii  (day  v  alues) . 

One  can  assume  that  the  particles  present  are  essentially  of  the  following 

three  species:  electrons,  ions  and  neutral  molecules  and  thus 

4 

m ^ m i  =  ^i8  x  10  . 

In  thiscase  (  ref.  2) 

ni  =  n^  =  10U  part/cm1  (day  values,  e,  electron  charge) 


Qi2  and  Qn  are  respectively  (Ref.  21  to  lb)  ; 
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volt/ cm 
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This  value  is  of  the  same  order  of  magnitude  as  the  corresponding 
See  beck  coefficient  for  semiconductors  (  Ref.  20)  whose  indicative  values 
are  reported  below: 

Material  Seebeck  coeff. 

Ge  (n-type)  853 

B  j  T(, 

h  c  3 

Si  -  Ge  -  alloy 
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3 .  -  The  electric  conductivity. 

The  expression  for  the  electric  conductivity  (  29)  can  be  simplified  under 

the  assumptions  of  the  Lorentz's  imperfect  gas. 

It  can  be  shown  that  the  following  relationships  between  the  coefficients 

L.'.'  and  the  transport  coefficients  hold: 
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nr 
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n,  m*  nu  Du 
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KT(ni  Dn  +  u,  0j2) 
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nj  n2  m,  m,  D, ,  I)2, 

1  L 

KT  (m  Dm  +  n,  Dl2 ) 

L" 

m  m/  D, , 

Pertinent  orders  of  magnitude  of  these  parameters  are  shown  in  Table  Ill. 

When  the  proper  substitutions  in  the  expression  for  0  are  made  one  finds 

that: 

,1  _  i  1  L  ^  '*°  ^oi  j  i  L>u  it  I'oo  ^ 


o  =  i!,t  -  '°  ,LV 

-e  tit  i 

I'oo  Ln-‘-.o 


1  1  1 1 
*'oi  ‘Ml 


I  1  i  1 
^Ol  ^O  i 


I  ' 

1,’,  .  -  --£>■-  =  0 
I’m 


since  the  first  two  terms  in  the  squared  bracket  are  of  order  one,  while  the 
third  term  is  of  order  (  6  °  )  .  The  vanishing  of  0  is  to  within  terms  of  order 
(  66  )  . 


This  important  result  indicates  that  in  the  present  case  the  plasma 


behaves  essentially  as  a  dielectric  medium  since,  in  the  steady  state,  a 
given  distribution  of  a  constant  electric  field  on  its  boundaries  does  not 
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originate  any  appreciable  current  but  produces  only  a  "polarization,  of  the 
medium,  as  indicated  by  the  induction  of  electron  (and  ion)  concentration 
gradients.  As  a  matter  of  fact,  it  is  just  the  presence  of  this  induced 
electron  concentration  gradient  which  causes  the  vanishing  of  the  electric 
current.  That  this  is  so  can  be  readily  realized  by  the  following  consider¬ 
ations.  As  shown  by  cc,s.(7)  the  electrical  conductivity  <J  for  an  isothermal. 

t 

uniform  three-fluid  plasma  would  be  simply  given  by  or,  for  a 
L.orentzian  gas,  by: 


which  is  obviously  the  Lorentz'  formula  when  the  ion  contribution  is 
neglected.  If  now  we  consider,  instead,  a  plasma  with  uniform  composi¬ 
tion  but  non -isothermal,  we  find,  according  to  eqs.  (7)  ,  that  the  conductivity 
(J,  would  be  given  by:  L.j2  -  (  'L'0(,)  .  But  the  second  term  is  of  order 

(  65  )  so  that  0,  is  practically  equal  to  Oj  .  Thus  the  presence  of  the 
temperature  gradient  would  not  alter  0  provided  the  electron  distribution 
remains  uniform.  If  however  we  consider  the  opposite  case  of  an  isothermal 
plasma  with  non-uniform  electron  concentration,  then,  in  the  steady  case, 
the  conductivity  0,  would  be  equal  to  Lj,  -  L^/L^.  That  is,  accounting 
for  the  relative  orders  of  magnitude,  it  would  vanish  identically.  Thus  it  is 
just  the  1  reaction"  of  the  electrons,  so  to  speak,  which  makes  the  total 
current  vanish  identically  in  a  first  order  stationary  state. 


CONCLUSIONS 
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We  state  the  more  important  conclusions  reached  in  this  chapter: 

0  For  a  ternary  macroscopically  neutral  plasma  in  thermal  equilibrium 
there  are  five  independent  state  parameters  but  only  three  gradients  of 
these  parameters  are  independent, 

»)  Chose  the  gradients  of  temperature  T,  of  electric  potential  6  and  of 

electron  concentration  C[.  as  such  a  set  of  independent  gradients.  Then, 
when  a  given  non-uniform  temperature  distribution  is  maintained  on  the 
boundaries  of  the  system  a  steady  state  is  eventually  reached  in  which  an 
electric  field  Vo  is  induced  throughout  the  system,  and  it  is  given  by: 

V<?  =  (l  V  T 

where  a  is  a  phenomenological  coefficient  depending,  in  general,  on  the 
locai  thermodynamic  state  of  the  medium.  Similar^,  an  electron  concen- 
'  ration  gradient  is  induced,  which  is  proportional  to  VT  according  to  a 
state  dependent  phenomenological  coefficient, 
in)  The  coefficient  q  ,  for  a  weakly  ionised  imperfect  Lorenlz-type  plasma, 
is  practically  a  function  only  of  the  degree  of  ionization,  increasing  with 
it.  alues  computed  for  an  indicative  case  (  ionospheric  F  layer)  show  that 
a  is  of  the  order  10  Volt/°K  for  a  degree  of  ionization  of  10'4. 
iv)  When  a  given  non-uniform  distribution  of  electric  potential  is  maintained 
on  the  boundaries  of  the  system  (  i.e.  when  the  plasma  is  imbedded  into  a 
tune  independent  electric  field)  a  first  order  stationary  state  is  reached  in 
which  the  electrical  conductivity  (7  will  depend  only  on  the  local  thermodynamic- 
state  of  the  medium  and  not  on  the  temperature  gradients,  say,  which  are 
induced  by  the  electric  field. 
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v)  For  a  weakly  ionized  imperfect  Lorentz-type  plasma,  CT  is  practically 
zero,  thus  indicating  that  the  plasma,  in  this  condition,  behaves  essentia  y 
as  a  dielectric, 

vi)  Points  ii)  and  iii)  give  qualitative  theoretical  support  to  some  experiments 


reported  in  the  literature.  (Ref.  11). 
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CM  APT  K!!  !V 

T11K  EXPERIMENTAL  1  IWESTIO  \']1Q N 

INTRODUCTION 

The  goner.* I  .iim  of  the  experimental  part  of  this  tad}  to  d*  tertiilno 
the  elect rodyn.tmu  properties  of  a  itommiforni  plasm.  In  particular,  the 
plasma  properties  of  air  about  a  super. '  ini*  vetii*  le  at  high  altitude  are  of 
in  to  rest . 

As  was  discussed  in  the  first  semi -annual  report,  this  investigation  would 
be  made  in  two  distinct  taps.  1  i  r  •  t ,  tie  experimental  apparatus  and  method 
would  be  tested  and  e\  ablated  using  an  >a:  .l>  generated  plasma  of  the  desired 
eh. i r.K  t e ri  ties.  O:  ly  it>  r  the  fir-t  - 1 •  •  p  1 1  t ^  been  satisfactorily  lompbted, 
would  the  mon  sop-ii-ti.  ted  mea-uren  ent  i  air  lie  attempted.  The  experi¬ 
mental  report  will  deal  primarily  witu  t  >  e.  il nation  of  tne  apji.iratus  and 
le  Unique.  \!  o  pn  l<  d  v.  ill  be  t  e  dat  ,  obtained  b\  utilizing  Inis  technique 
to  examine  the  pi  -m  ;,r  xlmtion  relax.it  >n  time  in  irgon.  Our  purpose  in 

presenting  t  i  data  on  argon  i  not  so  mu  n  for  its  intrinsic  value  (whi  h  is 
dubious  because  of  impu r it ie  ) ,  but  rat  ier  to  show  the  type  of  problem  and 
the  iiciur.ii  y  towniei  tne  technique  can  be  meaningfully  employed. 

The  experiment  1  t<  .  mi  pie  will  be  discussed  in  tile  following  order: 

First,  tlie  shock  line  will  be  considered  as  t  le  phi  sma  -  p  rodu  e  in  g  agent; 
second,  the  equipment  t->  me.i  re  t-  •  t  lei  t  rodynamic  character  ol  the  plasma 
will  be  discussed;  and  l:  .••!  ,  t  e  experiment  il  technique  wilt  be  evaluated 
using  the  data  obtained  from  experiments  run  using  an  argon  plasma.  This 
section  wi  LI  conclude  with  our  future  objectives. 
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t.  THE  PLASMA  GENERATOR 

Numerous  methods  of  generating  plasmas  are  available  to  tin  experi¬ 
mentalist.  Unfortunately,  at  this  state  of  the  art,  no  one  has  yet  devised 
a  plasma  generator  where  all  the  pr ope r tie.-,  are  controllable.  \t  bent,  it 
is  necessary  to  make  .»  compromise  between  tlu»  characteristics  important 

to  a  given  experiment  and  the  failings  of  the  plasma  production  system. 

« 

The  important  rerpii  ite  for  the  plasmas  in  this  experiment  is  that  the 
density  can  be  controlled  not  onh  in  magnitud<  but  also  spaeially.  Additional 
requirements  arc  Unit  the  gas  density  and  type  he  controlled  and  that  the 
plasma  can  be  contained  so  that  the  cle  tromagnetic  character  can  be  deter¬ 
mined.  To  these  ends,  a  shock  tube  as  a  plasma  generator  is  a  suitable 
compromise. 

It  was  shown  in  thi  fir  t  semi -annua  I  report  that  the  shock  tube  was  capable 
of  producing  known  plasma  densities  with  or  witnout  axial  nonuni  form!  t  ies. 
Further,  a  discussion  was  presented  as  to  the  i  ont reliability  of  the  plasma 
composition  and  lh<  u  ;<  of  the  shock  tube  .is  a  wave  guide  for  electrodynamic 
measurements.  Ns  is  usually  the  case,  it  i  not  the  general  behavior  of  .1 
device  that  is  important,  but  rathe  r  the  spe  itic  departures  that  must  be 
analyzed. 

In  a  conventional  shm  k  tub  ,  .1  pressure  ratio  is  establishes!  across  a 
diaphragm.  When  the  oiapnr.igm  either  bursts  or  is  broken,  the  thermo¬ 
dynamic  character  of  the  lug  of  air  behind  the  shock  that  is  formed  can  be 
computed  a  priori.  In  principle,  by  adjusting  the  initial  conditions  (in  both 
the  driver  and  driven  sections),  a  plasma  can  be  generated  with  specific 
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properties.  However,  energy  dissipation  processes  which  cannot  be  accurately 
considered  in  the  predicted  shock  calculation  may  substantially  attenuate  the 
theoretical  shock  velocity  in  the  driven  section.  \\  ith  a  decrease  in  velocity 
of  about  20%,  the  predicted  electron  density  in  the  plasma  may  be  reduced  by 
more  tlian  an  order  of  magnitude. 

To  compound  this  prediction  error,  in  a  shock  tube  where  the  diaphragm 
i  self-rupturing,  there  is  an  additional  un»  ortainty.  The  variation  of  mechani¬ 
cal  properties  of  the  diaphragm  may  introduce  an  additional  order  ol  magnitude 
of  uncertainty  in  the  cU  .  iron  density  value. 

Since  it  is  desirable  in  this  experiment  to  out.  in  plasma  densities  within 
a  small  range  for  e.nh  operition  of  the  >  tock  tu  n-,  the  following  techniques 
have  been  employed.  Fir  t  noth  tne  initial  driver  and  driven  pressures  are 
measured  by  mechanical  gauges.  A  Wallace  ami  Tiernan  0  -  50  mm  gauge  is 
Used  to  measure  tin-  initial  driven  presstin  .  Ft.  ortunatc ly ,  this  gauge  when 
used  below  2  mm  llg  requires  a  long  ti  h  bet  on-  an  equilibrium  value  is  ob¬ 
tained.  in  principb  ,  iiouew-r,  the  driver-driven  pressure  ratio  can  be  main¬ 
tained  to  within  3%  for  a  driven  pressure  of  t  mm  llg.  \itcr  calibration,  tue 
absolute  value  of  the  initial  driven  pressure  is  al.  o  accurate  to  about  3%  at 
1  mm  llg,  improving  to  .  3%  at  I  )  mm  llg. 

It  might  be  noted  Una  many  shock  tulles  am  operated  using  metallic 
diaphragms  where  the  rupture  pressures  arc  predictable  to  within  5%.  However, 
it  was  deemed  in  pi  r.ttivi  lli.il  no  met  illn  fragments  cxi  t  in  the  shock  tube 
since  these-  might  damage  the  tube  and  mil  rowave  equipment  and  possibly 
produce  erroneous  elect  rodynumic  measurements. 
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In  order  to  m. tint. tin  the  driver-driven  pressure  r.ttio  precision,  it  is 
necossa r\  to  break  rather  than  pressure  rupture  the  main  diaphragm.  In 
tue  lirst  semi  annual  report,  a  method  ti-ini*  hot  wires  stretched  across  l!ie 
diaph raym  was  desc  ribed.  IVhih  mi  method  produced  satislae lory  breaks 
wit  i  the  conditions  Hie  tube  was  then  ope  ratine,  tue  method  was  cumbersome 
since  tue  wires  recpiired  replacement  alter  each  -hot.  In  addition,  teaks 
developed  in  the  <  ompli  i  ated  uomnetalli  diaphragm  holder.  Further,  since 
in  the  future  it  is  t  \pe«  ted  t.uat  t  »c  ten  pi  r.itnre  ol  me  driven  tube  would  bt 
elevated,  to  Is  were  run  t  »  lind  .»  pi  «  sti  -..lit  » h  l  <  it  h  i  _;h  teinperatu  res  that 
could  he  broken  u-ii  c  to.-  >t  wir  le  .  \'o  such  plastic  could  h«  found. 

1  he  hot  wire  technique  was  abandoned  in  I  c  cor  cl  a  mech.ude.il  *  tilting  system 
descrined  below. 


a  Nil  roai  n  .  ii  mi  be  r 
I)  Sec  Olid  diaphragm 
c  Piston 
d  Stops 


e  n  r i v •  r  tube 
1  I  tod 
j  Cutter 

M.i i  it  diaph ra  i>m 
i  j)  riven  tube 


Figure  I.  Diaphragm  breaker 
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Thu  main  diaphragm  consists  of  four  slu-cls  of  .  010-inch  I.exan*.  This 
combination  has  a  rupture  st  reugtii  of  about  JOL)  psi  for  the  three-inch  tube 
diameter  in  this  experiment,  til  normal  operation,  the  pressure  differential 
across  tiie  diapnragm  is  about  150  psi.  As  shown  in  Figure' l  above,  a  fonr- 
bladed  pointed  cutti  r  i>  attached  to  a  quarter  inch  rod  that  extends  the  length 
of  the  driver  tube.  On  t  o  oppo  ite  end  of  the  rod  from  the  cutter,  a  piston 
is  attached,  i'he  pi  t  »n  opt  rates  in  .«  chamber  which  is  divided  by  a  second 
diaphragm  consisting  of  four  sheets  of  l.ev  m.  1'n  break  the  main  diaphragm, 
the  section  of  the  chamber  not  inc  1  tdhi„  tin  pi- ton  is  pressurised  with  nitrogen 
until  the  second  diapnragm  rupture  .  Pressure  is  then  exerted  on  the  piston 
which  a<  (  elel'ate  t  ie  c.Utter  attaining  a  Vel  »  its  sufficient  to  CVlt  the  main 
diaphragm  in  about  one  millisecond. 

This  system  has  proved  highly  successful  but  slightly  cumbersome  because 
of  the  Second  diaphragm.  However,  it  offers  tne  advantage  that  any  reasonable 
driven  pressure  can  be  list’d  since  the  tnickrwss  of  the  main  diaphragm  is  not 
ritical.  Further,  it  allows  t  .<.  u-e  of  any  plastii  sheet  that  can  in  cut,  thus 
permitting  the  use  of  high  temperature  material-*.  Modifications  of  this  system 
are  now  in  progress  so  as  to  replace  the  second  diaphragm  witii  a  monostable 
valv  e  assembly. 

The  present  shock  tu  k  y  tern  operating  at  an  initi  >1  driven  pressure  of 
about  l  mm  Hg  and  driver  pressure  of  150  psi  lias  an  average  deviation  from 
the  mean  running  conditions  of  about  70  tt/see  at  7000  it/sec  shock  velocity 
or  about  1%  mean  deviation  in  velocity. 


Produced  by  General  Fk  trie  Co. 
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In  argon  at  7001)  lt/sec,  this  deviation  in  vi  Unity  corresponds  to  about 
30 n/a  in  the  equilibrium  electron  tlensity.  Ii  should  be  noted  th.it  these  figures 
are  obtained  from  the  exp.  rimont.il  running  conditions  of  this  sho  k  lube  and 
represent  the  menu  shot  to  shot  variation  that  »  an  be  expected.  If  the  theore¬ 
tical  equilibrium  electron  density  calculations  are  made  using  the  actual  shock 
velocity,  the  initial  dri.en  pressure  and  tne  gas  composition,  tin  equilibrium 
electron  density  of  any  given  shot  »  an  he  determined  to  about  3fo.  This  i 
because,  for  eacn  shot  t  c  vekuity  ( wl.i.  a  is  lie  important  vari.tbk  )  can  be 
tlete  rmiued  to  about  .  I  ’  j  using  tne  relit  ted  n.i  ruw.ivi'  sign.  1.  M  i  <  well  known 
tit  tt  in  the  particular  ease  oi  argon,  the  approach  to  equilibrium  ioni/i.ition  is 
-low.  Thus  t  to  equilibrium  calculation-  nm-t  '*•  used  only  as  an  indication 
ot  the  asymptote  vain  ot  .  leitron  den-itv  behind  the  shock  it. I  e  iTu  c  e . 

It  lias  been  shown  l  at.-  far  t  tat  the  shock  lube  is  expefime.n tally  capable 
c>f  generating  known  a.,  rage  plasma  den  ili  -  in  a  repeatable  manner . 

Another  source  of  •  rrors  in  the  ini.  rpret  lion  of  the  experimental  results 
is  the  nonu.Tfortn  radial  ili  trituitiuu  oft  te  plism.i  properties  dm  to  boundary 
layer  effects.  For  in.-tuae,  tne  estimated  raui.d  nonunifo rmity  for  air  is 
shown  in  Figures  J.  anti  ).  in  l  e  present  series  ol  experiments,  the  microwaves 
propagate  behind  l  a  -hoik  into  rfuce  at  a  di.  lance  where  thi  distortion  of  the 
radial  distribution  ot  thi  <  l .  ct  romagneti .  ti>  Id  l  till  negligible.  1  ill  t  lese 
effects  were  not  taken  it  to  a.  count  in  tne  interpretation  ol  the  experinu  nt  il 
results.  At  the  present  time,  .  alculatioiis  are  being  .ondmted  to  include  the 
radial  nonuniform  di  lid  ml  ion  in  the  analyse  ot  the  wave  guide  fundamental 
mode. 

The  pre  ent  exp«  riieent.tl  iuve  ligation  is  c  outined  to  the  analysis  of  the 
axial  I  y  non  uni  I  on  n  <  I  -  I  em  density  in  argon. 
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The  controllability  of  the  density  gradient  in  air,  as  described  in  the 
first  report,  lias  not  yet  been  utt<  mpted. 

“•  El'ECTKOM  AGNKTIC .  MKASUKK M K NT S 

In  the  previous  section  the  sliock  tube  was  described  as  a  plasma  generator. 
The  salient  feature  of  this  experiment  is  that  t  e  shock  tube  also  can  be  used 
as  a  wave  guide  whose  properties  depend  upon  the  electromagnetic  properties 
of  the  generated  plasma.  1  lie  range  ol  frequent  ies  that  can  be  propagated  in 
the  fuudainei  t  il  transmission  mode  1  hl(  is  small,  which  limits  the  lange  of 
electron  densities  tii.it  e.tii  be  measured.  llowevi*r,  it  will  be  shov.ii  th.it  \iithin 
this  range,  a  remarkable  de  ree  of  prn  i.  ion  can  be  obtained. 

It  is  will  known  that  a  cylindrical  tube  can  ait  a  a  transmission  line  tor 
certain  frequencies  ot  electromagnetic  waves.  T1  t  le  medium  in  t.  e  \v  is  e  guide 
is  a  vacuum,  for  tin  1  l'.lt  mode  t  le  cut-oil  frequeney'  is  given  >y 

l.  8  1  i 

t:  ”  'fit 

where  K  is  the  tube  radius  and  i  i  t  ie  vt  locitv  of  light  in  a  v.ieuum.  \ll 
frequencies  abovi  this  "cut-oil"  .trequemy  -  m  hi  propagated,  but  the  geometrii 
field  configurations  In  come  very  compii  ale  .  Since  it  is  necessary  to  measure 
thi  electric  and  magneti  iields  in  tiiis  cxpi  rimeiit,  it  is  desirable  to  have  a 
situation  where  only  one  mode  can  propagate  in  tiie  wave  guide.  I  bus,  only 
frequencies  between  t>  e  ■  ut-off  froquen  y  given  >y  the  above  equation  and  the 
cut-ofl  t  re  quern  y  for  the  next  modi  an  used.  This  upper  frequency  limit  then 
is  given  by 

-ID5  c 
"  ;<» 

which  corresponds  to  the  cut-oft  frequency  of  the  TMo1  mode. 
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For  our  shook  tube  with  .1  radius  of  C  75  cm,  the  t  utoff  values 
establish  the  frequency  r.iu^c  from  i350  to  3Uo5  mc/scc.  Furt n  r  limitations 
on  this  range  are  impos.-d  by  the  inability  to  turn  the  system  as  the  low  fre¬ 
quency  limit  is  approached  because  t  ie  w a\  t  guide  wtvc  length  becomes  Very 
long,  extending  to  infinity  it  cutoff.  Fi  itr<  I  shows  tue  relation  between  the 
wave  guide  waveleiiglu  tnd  tue  generated  trequ*  ncy  tor  this  shock  tube.  The 
propagation  limits  tor  tin  I  K,|  modi  are  siiow.i.  Whei  tin  wave  guich  wave¬ 
length  exceeds  about  I  )  >111,  tin  system  bci  ome  s  lifli  ult  to  tune.  lh  me  the 
practical  operating  r.m^i  is  about  i-loU  toSd  im  /sec. 

\\  ion  the  w.»v<  guide  or  shock  tube  •  out  tins  a  uniform  plasma  the  propa¬ 
gation  depends  upon  t  >  plasma  1  reqiie’  cv  and  the  electron  collision  f roquency 
in  lh«  plasma.  in  tue  ideal  case  vv  iftre  t  a  crdlt  ion  1  roquency  is  assumed  Jtero, 
t  ie  propagation  1  on 1 , . ; .  t  |  ^  iui'  t.ic  i  ,j|  11  od*  is  given  ia ; 


a  i  ■  re  j  is  t  tie  a  n  gu  1 .1  r  t  re  pie  lie  y  ot  t  .  >•  1 1  .  t  r  01  n  a  gm  •  1 1 1  wave  and 

“O  j  j  '  •'  t  •'*'  plasma  1  r  1  |  iii  in  y  .  is  ri  lit'  d  to  t  :u  •  ilci  tl*i>li  del  ~  it  y  u, , 

by  tile  equation; 

^p  *  \  U  *  /nl%  A  (i) 

where  n^,  mi  «isur<  llu*  number  ot  i  loi  t  ron.i  per  cubic  cuntimcti  r.  From 
equation  (I),  l  ie  cutoff  frequency  u>c  is  given  by 

J  L 0 

-J,  -  ..  .  I  7 1  x  Id  =  u,- 

P 


(3) 


iwmiiawj 


•;•,'*■ 'TW'fn' 


[UiiU: 


ilHiiiiiiiiffili 


iirajiiiiliilA1 


S‘iHH 


qun; 


mm  sik 


MimHHR 


■mnu 


t 


Frtqutncy 
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u)  . 


3.  18  x  10 


9 


U) 

-  -  o.  83  x  10  =  no 


(9) 


u)t  aiUbt  be  sm.illi-1'  than  tin:  cutoff  angular  frequency  of  the  TMe|  mode  or 
the  wave  guide  ahead  o'  the  plasma  inti  flat  e. 

Therefor,  ,  v.itl  t  i,  maximum  in  qwency  5  >3  >  nic/sr,  ,  the  upp,  r  limit 
of  electron  density  at  whi*  the  c  le<  t  romajjm  tii  Held  call  propagate-  in  t  ie 


plasma  is: 


Id 


n  3.3x10  e  le  l  t  rail-  a  l 

I'niaX 


(3) 


ID  mea  su  n  tie  pr  .perties  ol  tae  pl.sma,  two  types  of  tii  Id  detectors 
are  used.  I  ir  t,  t  ,  el.  trie  fi.  Id  d,  t  tor  is  n  antenna  wniui  i-  mounted 
flush  with  the  w  .11  H  t  a  snoek  tone.  1  i  ant.ona  is  sonneted  to  a  mi.  ro- 
wave  diode  type  IN  33.  H  the  detector  op,  rating  in  the  "sMn.,r,  law"  reuion, 
tv  output  volt.,*,  tram  t  ..  mod.  is  pr  -p-rtm,  l  to  the  square  of  the  ,v.-ra*e 
ele.trii  field  at  l  ■  ,11.  Pi,  maximum  «>  tput  from  a  typical  d.  toe  t  or  for 
an  input  power  >  t  3  )  milliwatt-  is  about  1  >  n  illicolts 

The  nia.net.,  d.  t.  tor  i-  like  me  el.  tri,  li.dd  dete,  tor  except  that 
h.  tween  the  vail  of  t  tub.  and  mtenna.  loop  is  mad,  .  The  pl.n  e  of  th, 
loop  lies  in  apian,  inc  l  .dm*  .  „•  axis  Of  the  tube  ,t,d  the  microwave  antenna. 

A„  ininp.i  . .  a.,, . Ur,.  -Mr., or  i-  . . .  «  ^rno.r  ol 

tn,'  loop  is  about  .9  ,u. 

\  spurious  siend  was  noted  on  the  <  Petri  probe  early  it  the  investi¬ 
gation  (shown  by  the  arrow  in  Fiuure  5  ).  Subsequent  investigation  without 
the  microwave  signal  ,  learly  showed  that  a  current  was  generated  in  the 
.  ry8Ul  probe  w  i  h  was  directly  associated  with  the  shock  arrival  (see  Fig.  6). 
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Figure  5. 


Shot  Number  <£3 

Sweep  Speed  50|isec/cm. 

Wavelength  ZZ  cm. 

Upper  Trace:  Pressure 
Detector 

Lower  Trace:  Signal  from 
microwave 
detector 

Note:  White  arrow  indicates 
interference  signal. 


Figure  6. 

Shot  Number  Z 4 
Sweep  Speed  50  (isec/cm. 
Upper  Trace:  Pressure 
Detector 

Lower  Trace:  Interference 
Signal. 


Figure  7. 

Shot  Number  51 

Sweep  Speed  100  p,sec/cm. 

Wavelength  ZZ  cm. 

Upper  Trace:  Signal  from 
electrostatic  probe 
(Amplification  l  V/ cm) 
Lower  Trace:  Signal  from 

Microwave  detector 
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A  hii>h  impedanci  probe  w.«s  made  and  insi  rti  d  into  the  shot  k  tube  m  pi  te  e 
ot  the  crystal  detector.  Following  tests  with  ind  without  microwave  signal, 
typical  voltage  curvi  -  like  th.it  in  Figure  7  v.crr  obtained. 

Cali  illations  conducted  for  a  ste.itly  tati  condition  of  a  uniform  plasma 
in  contact  with  i  nut,. Hi,  w.,11  indicate-  th  t  a  potential  differetu  r  o0  i;iviii  by 


K  1 

t  • 


(o) 


mist  be  established  in  or  hr  to  have  a  zero  eb  Iron  i  nr  rent  from  the  plasma 
to  the  wall.  With  the  valiu  of  i  on  -  elei  t  r  >n  mas-  ratio  (  ^  *  *N  for  nvon  and 

M  ~J 
•  ' 

an  electron  temp.  r.  tin  e  (i;)  on  the  order  .f  1  )*  °K,  equation  (-)  .ies  a  value 
4*o  19'’  faiiite  of  i  Volt  .  I  hi  ~  vain,  mini  ides  with  the  measured  potential 

on  tlie  eh  tro:  tati.  probi  in  l*  i  ’un  <  .  I  d  ii.ti  rmrcii  e  sianal  was  removed 
oy  introdu  in c  a  small  p.irillil  indu  tain  •  on  tin  RF  siii  of  t  u  el.  trie  field 
probi  . 


9  ae  c  otipltii t;  of  tin*  n  icruwave  strnrr  iloi*  to  the  w  ive  ^uide  is  an  import. ml 
consideration.  Tim  spe,  i ! i .  coupling  condition  i-  that  ii  the  wave  t-uide  was  of 
infinite  length,  tne  rcflecteo  pi > v.  c  r  1 1 1  r i > i i li i ■  tie  a1  tenua  would  be  /. <  ro.  Stated 
dilfe  l  vntly,  il  a  in  e  t  a  1 1  i  pi  ton  is  pl.t  c  ed  ii  tin  \\  .  v  i  ■  ^  w  i  1  e ,  a  .  t  inline  wave 
pattern  would  be  establi.-  led  such  that  linear  lispla  ement  ol  the  piston  would 
result  in  a  proportional  linen-  si  ift  in  a  mil  i  an  of  the  standing  wave  pattern. 
Or,  if  an  electri  field  1<  te.  tor  i,  statioiiai-)  to  the  wave  t>uide  and  the  piston 
is  linearly  displaced,  t  u-  maumtude  of  the  I'i .  I  1  will  vary  sinusoidally  with  the 


piston  motion. 


9.d 


In  the  previous  section,  Hie  ciiar.ti.li  r  of  the  plasma  generated  in  tlie 
shock  tube  was  described.  For  the  amplest  i  ppro.xinia  tion,  the  generated 
plasma  could  be  considered  as  a  uniforn  ly  ionized  slug  of  gas  radially  filling 
tin  tube  and  extending  about  oik  meter  behind  the  shock  front.  The  electron 
density  of  this  plasma  detint  d  by  nt.  con  I  >.•  t  .  ken  as  the  tiieri'  odyiiainit 
equilibrium  value.  If  nt.  i-  so  large  that  the  propagation  constant  k,  in 
equation  (1)  bi  comes  .,n  imaginary  qu.u  tity,  tile  shock  front  acts  like  a  reflect¬ 
ing  piston  The  signal  from  a  stationary  electric  field  detector  in  front  of  the 
shock  wool  1  vary  luusoid.ctly  with  time  as  !:n  shock  moved  towards  it. 

The  time  required  to  i  omph-te  one  sine  cycle  in  tin*  detector  signal  is 

related  to  the  shock  \el>».  it y  by 

\u 

'  ~~2~~ 

where  \  is  tm  wave  guide  waved  ii.tli  .ml  v  is  i  e  shoe  k  veloc  ity . 

Since  \  ,  and  7  c.ii.  easily  be  dutern  tied  to  ijett.  r  than  0,  1%,  the  shook 
veloc  ity  eall  also  measure  to  at  least  o.  If,. 

In  f  igure  Ha,  a  vrl  u  ity  history*  (aii  iie  determined  a  the  shock  moves 
toward  t  ie  detector  it  t  e  end  of  the  tu  «  .  T  .is  history  is  showi  in  Figure  8b 
and  is  typical  ol  sho  k  line  b  iiavior.  in  1  igure  Ha,  toe  m  rease  in  amplitude 
ol  tile  detected  signal  is  din  to  Hit:  fact  tli.it  the  i  lierowave  signal  is  attenuated 
in  the  wave  guide.  Hence,  t  ie  standing  wave  ratio  decreases  as  the  shock 
interface  distance  from  ti  c  detector  increases. 

As  the  shock  passes  ttie  detector,  the  character  of  tin’  detected  signal 
changes.  To  illuslr  te  t  sc  changes,  t  o  set  of  graphs  in  Figure  9  have  been 
prepared.  In  all  cases,  the  time  average  of  the  transverse  electric  and  magnetic 
field  at  the  wave  guide  surface  are  shown.  The  dotted  Lines  represent  field 
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conditions  of  the  previous  cycle  so  as  to  contrast  the  change.  If  the  reflecting 
surfaces  (piston  or  plasma)  are  assumed  moving,  the  field  conditions  arc- 
drawn  in  the  frame  of  reference  of  the  moving  system.  Hence  a  stationary 
detector  would  measure  the  field  distributions  in  the  reverse  order.  Since 
tlu-  experimental  data  hive  all  been  obtained  using  argon  shocks,  the  plasma 
examples  in  Figure  9  illustrate-  field  conditions  for  various  .axial  electron 
density  gradients,  assuming  some  absorption  inside  the  plasma. 

In  the  case  (a)  of  Fig.  9  the  penetration  inside  the  metallic  piston  is 
very  small  diminishing  to  zero  as  the  resistivity  approaches  zero.  Note  that 
the  transition  in  the  transverse  magnetic  field  occurs  at  its  maximum  and  the 
converse  occurs  tor  the  electric  field.  Case  (  d)  for  a  steep  plasma  gradient 
is  similar  to  (a)  except  that  the  penetration  depth  is  greater.  A  typical 
experimental  example  ot  this  i  ,isr  is  shown  in  Figure  10.  Similarly  with 
case  (gl,  where  the  plasma  gradient  is  less  steep  than  case  (  d)  tin-  penetra¬ 
tion  depth  has  increased  and  tin-  departure  points  on  the  electric  and  magnetic 
field  increase  and  decrease  respectively.  Note  that  the  electric  field  inside  the 
plaso  a  decreases  with  a  rate  changing  exponential  since  the  electron  density  is 
increasing,  experimental  behav.  ior  of  this  type  is  shown  in  Figure  11.  An  en¬ 
largement  ot  Figure*  11  has  been  made  in  Figure  1  i  where  the  detected  electric 
field  inside  the  plasma  and  the-  clectri<  field  of  the  previous  cycle  are  compared. 

Figures  1  1  and  H  were  taken  so  that  it  would  be  possible  to  discern  the 
departure  points.  These  pictures  were  made  using  two  electric  field  detectors. 
One  amplified  detector  signal  was  fed  to  the  horizontal  deflection  plates  of  a 
cathode  ray  tube,  .he  other  to  the  vertical  plates.  Spacially,  these  crystals  were 
located  several  wavelengths  apart  on  the  shock  tube.  The  trace  becomes  un¬ 
blanked  about  three  cycles  before  the  shock  passes  the  first  detector.  Hence, 


l 


i: 

i. 

i 


i 


Figure  10. 

Shot  Number  179 

Sweep  Speed  20  psec/cm. 

Wavelength  18.9  cm. 

Upper  T  race:  Signal  from 
electric  microwave 
detector 

Lower  Trace:  Signal  from 
magnetic  microwave 
detector 


Figure  1 1 . 

Shot  Number  116B 

Sweep  Speed  50  psec/cm. 

Wavelength  44.4  cm. 

Upper  Trace:  Signal  from 
electric  microwave 
detector 

Lower  Trace:  Signal  from 
electrostatic  probe. 
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Figure  13. 

Shot  Number  110 
Wavelength  44.4  cm. 
Horizontal:  Signal  from 

microwave  detector 
Vertical:  Reference  signal 


Figure  14. 

Shot  Number  1 1  8 
Wavelength  44.4  cm. 

Horizontal:  Signal  from 

microwave  detector 
Vertical:  Reference  signal 
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the  normal  standing  wave  sinusoid  is  established  is  an  ellipse.  (The  non¬ 
symmetry  is  due  to  the  nonline  rity  of  the  detectors.)  As  the  shock  passes 
the  first  detector,  the  departure  is  clearly  evident.  The  location  of  the 
departure  relative  to  the  shock  front  can  he  made  by  using  the  corresponding 
linear  time  scale  pit  lure. 


Figure  15. 


Referring  to  Figuri  15a,  note  that  if  tae  departure  occurs  on  the  outside 
of  quadrant  L  or  3,  the  function  is  decreasing  faster  than  the  sine  function  and 
presumably  denotes  the  rate-increasing  ty-pe  of  exponential  function.  However, 
if  the  departure  is  oppo  iti  as  shown  in  Figure  15 b,  the  detected  signal  in  the 
plasma  is  decreasing  slower  tnan  the  »ine  func  tion  and  corresponds  to  the 
case  described  below.  The  arrows  denote  the  direction  of  trace. 
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Case  (j)  of  Figure  <»  represents  a  possible  field  configuration  if  a  slow 
axxal  gradient  exists  behind  the  shock.  In  this  case,  the  critical  c  lectron 
dc-nsity,  where  the  propagation  constant,  k..  is  ,.ero.  is  located  a  long 
distant  e  behind  the  physical  sho,  k  front.  Hence  there  is  a  region  of  the 
plasma  behind  the  shock  interface  wnen  : 

r - t0 

up  v  **  -  --17  *  10 
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of  tlu-  M«„„l  whirl,  appro., rim,,  «•  ro  i„  ,!„■  n  rion  ,„  hi«h  rl,,,™ 
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Figure  16. 


Shot  Number 
Sweep  Speed 
Wavelength 
Upper  Trace: 
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50  |iset/cm. 
44.  4  cm. 
Signal  from 


microwave  detector 
Lower  Trace:  Signal  from 
electrostatic  probe 


Figure  18. 

Shot  Number  164 

Sweep  Speed  20  p.sec/cm. 

Wavelength  22.49  cm. 

Upper  Trace:  Signal  from 

electric  microwave  detector 
Lower  Trace:  Signal  from 

magnetic  microwave  detector 
Blanking:  Pressure  detector 
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equation  (  i)  assuming  the  collision  frequent y  zero,  lienee  for  each  micro- 
wave  frequency  used,  t  electron  density  at  the  ’’reflecting  poii  t  '  is  known. 
Thus,  in  principle,  hy  repeating  the  sairu  slio.  k  and  changing  th «  microwave 
frequem  y ,  .1  plot  of  tit  a  ran  den  ity  vs  distance  behind  the  shock  can  be  made 
lor  the  range  of  dcr  it  it-—  described  in  Section  Z. 

This  general  prat  i  dt.ro  was  used  to  e  t  t  ij  li  sh  the  density  p  r  of  i  I  •  behind 
an  argon  shot  k  for  sec  era  I  shock  v<  lot  ties.  However,  rather  than  change  the 
f r<  plenty  between  eat!'  snot  (w  i  it  is  a  laborious  protedure  since  the  system 

must  be  rt  tuned  ,  .1  series  of  test  w<  re  omhuted  at  ea  h  frequency. 

To  mark  1  u-  p  y~ic«il  arric.it  <>1  tlit  shock  a  1  ry  t  tl  presstir  probe  was 

mounted  flush  with  t  hi  tube's  surf  ice.  1  >•  output  from  this  crystal  was  fed  to 

a  pulse  criiuit  whit  h  in  turn  1  nt  of!  the  Wear,  ii.lt  ns  ity  of  the  measuring  scope 
for  5  mi  rose contis.  \  type  at  trace  of  tin  t  lettrit  and  magnetic  probes  is 
shoevn  in  Figure  18. 

Note  that  a  di  s  pit  t  t-niei  1  or  re*  t  ion  mu  t  be  made  betccet  n  t  :io  signals 
since  tin  probes  are  lot  a  ted  3  m  apart  along  the  axis  of  the  tube.  This 
correction  is  mailt  by  lir  t  determining  t.a  vein*  ity  as  first  ribed  in  Set  lion  Z. 
Pdnce  tae  distinct  betevt  en  proses  is  know;  , by*  using  the*  measured  velocity,  a 
time  torn  ction  of  events  i<  made,  and  tn.  di.-l.iiue  beteveen  shock  arrival  and 
''reflecting  point ''  tail  In  determined. 

Using  tn i s  procedure,  Figures  19,  ZD,  and  Zl  were  obtained.  From  these 
curves,  a  eomposit  curve  can  be  drawn  showing  the  electron  density  as  a  function 
of  distance  behind  the  shock  for  various  shock  velocities  (Figure  ZZ)  .  The  de¬ 
gree  of  randomness  in  the  experimental  values  is  due  to  the  uncertainties  which 
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wore  previously  discussed.  Acta. illy,  t  h  <  ■  result  s  are  surprisingly  good 
taking  into  account  that  t .  i  *  ionization  rate*  of  argon  is  a  v<  ry  sensitive  function 
of  impurity  content. 

S.  i  the  purpose  ot  t  a  i  s  particular  investigation  was  to  demonstrate 
tile  capability  of  the  measuring  toclmiepic,  no  special  precautions  regarding 
tie  impurity  ot  tile  argon  v. a  re  taken.  It  might  ie  noted  til  it  t'  i  data  presented 
was  obtained  by  using  argon  from  t  ie  same  tank.  Upon  (hanging  tanks,  a  com¬ 
pletely  elide rent  set  ot  >ur.es  v. as  obtained.  I  he  estimated  impurity  content 
-  i  i 

is  between  Id  and  Id  tor  tue  data  pre  sent  e  t . 

The  theoretic  a  I  i  p.i  I  i  nrium  electron  densities  for  argon  arc’  computed 
as  an  appendix  tat  is  lion.  From  t  ie  •  results,  asymptotes  un  lie  located 
on  the  curve  s  in  It  ires  1  »,  2J,  I,  and  J.i  '1  ie  relation  shown  between 
measured  and  t  leoreti  I  .  h  ctr  on  dem  ities  indicates  the  hi  gli  precision  which 
cat.  tie  attained  wita  t  i  ~  l  •cntltpie. 

An  inde  pende  at  correlation  ot  tac  m  r »  salts  ts  presetited  in  P  igurc’  «  i . 

T  he  re  I  exation  time  defined  by  tin  other  it  a  '  tig.itors  refers  to  the  degr<  e  of 
ionix.it  i  oil  as-  me  a  s  u  r  •  a  1  by  t .  1  It  i  no s  .  i  v  of  t  n  i  plasma.  F  h  ct  s  e  a  u  t  tots  claim 
a  direct  relation  between  tac  I  >  miinos  i  t  ■.  and  the  electron  density  .  J  he  point 
plotted  from  our  data  repn  sent  a  relax aliot  time  Ic'fined  as  t  ie  time  tequiled 
for  the  plasma  to  attain  ot)%  of  its  ecpii  I  i  ar  i  am  \alue. 

Bec  ause  of  tiie  possiale  ambiguity  in  the  detinilions  of  T,  only  an  ordc  r 
ol  magnitude  agreement  t.  nit  tiling  fill,  I  the  results  tor  itn  pu  i  ity  c  on  t  e  n  t 
in  the  range  of  M  3  to  I  )  1  are  extrapolated,  such  an  agreement  is  clearly 
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•1.  FUTURE  WOKK 

Since  tile  ar^on  shock  lias  boon  den  on  t  rate  1  to  be  i  simple  and 
predictable  method  of  obtaining  ail  axially  nonuniforni  plasma, flirt  icr  ex¬ 
periments  will  at  made  using  argon  to  t  in  t  k  the  analysis  of  Chapter  V. 
When  satisfactory  correlation  of  this  analysis  has  been  accomplished,  our 
investigation  will  be  dire  ted  towari  plasma  e  radiants  in  air  shocks. 


APPKNLF'  A  TO  CIIAPTLTv  IV 


r  *  *  IMP 

I  1  I 


KQmLimm_!.\i ohock  c  \lcul ation  for  pukf  MiGON 

Equilibrium  conditions  behind  normal  shoi  k  waves  in  argon  have  been 
computed  using  a  nutuxl  o  it  lined  i»y  . .  1 , .  R<  slo  r ,  S .  C  ,  Lin,  and  \  .  Kant  row  i  t  /, 
^  *1  hfe  Production  of  Iiiyh  1  einpe ratu re  Gases  in  a.,ock  Tubes,  Journal  of  Ajtplitd 
P  vsi.s,  Vol.  .f  i ,  \o.  1J  (Do  ember  I  * 5 -  1 '  ie  i  ,il  illation  consiib  r-  sinyle 

>iii /.at ion  of  t  he  a  r  got  at  c  o,  but  e  lei  t  rot, i  exit  t  i>>n  and  trai  It i  pie  i  oni /,<i  tion 
.ire  ignored.  The  assini  ,  ion-  u-ed  int  radio  <  ,  rror.-  of  I  L'%  at  15,000°K, 
and  are  i  on  side  r.i  nl\  mi  r.  i‘iur,ite  at  lower  temperatures. 

I  In-  conservation  equations  for  mass,  momentum,  and  ei  rgy,  plus 
t  ie  equation  of  state  aid  pertinent  the  rmodytiumii  relationships  ha',  e  been 
written  tor  aiqon  undcr^uiin;  tie  ioni/.ation  reution; 

\  l  e 

I  ae  sunsi  ripts  1  ,.n<i  J  denote  condition  in  front  of  and  behind  t  ie  shock,  and 
the  in  it-  used  a  re: 

density,  p,  -lugs/ft* 

pressure,  p,  atm 

v<  loeitv,  u,  ,t ,  ft/ sei 

t  ernpe  r  it  u  n  ,  1' ,  °K 

enthalpy,  h.  It  "/sec" 

elec  1  ron  den-ity,  iq  ,  parti,  les/.  in' 


and  a  is  tiie  traction  il  degree  of  ioni/.»lion  of  ,ir:on. 


Pi  Mg  =  Pi  (Ms 


u 


1  1 


“i 


> 


•3116.8  p,  +  p,  af  M“  =^116.8  p,  l  p2  aj  (Ms  - 


hl  !  7  -M  M.;  =»i.  »  -1  a“  ( M  s  -  —  —  ) 

.i, 


3  u 

h  -  5.637  x  10  ( l  t  (i  )  'I  f  4.075  x  l J  a 


a  - 


-.510  x  1 J5  _!>  exp  (±*=  iii  -‘i’5- 
T5'“  \  T 


1  -  I 


7  - 


p  =  l .  0t>5  ( 1  •  it )  p  i 


a,  =  6t.  j  jy- 


n  1 1 

n  .  =  7.  3  1  x  I  J" - JL 

1  t  a  T 


An  iterative  e  lution  scheme  was  pro^r  mimed  for  the  Bendix  G-  15 
diait.it  computer  u.-im>  the  ove  r  relaxation  method.  Thus,  given  input  con- 

f  . 

ditions  Ms>pj  ,  and  1  j  ,  i  unditions  behind  the  siiock  nny  In-  obtained. 

Figures  .’.-1  and  .15  present  some  cat  t  ted  e lei  Iron  densities  as 
functions  of  shock  velocity  and  initial  pressure  .  file  initial  temperature 
was  300°K,  and  it  was  tound  taat  tin  same  results  are  also  Valid  for  some¬ 
what  tower  initial  temperatures. 
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CHAPTER  V 


ANALYSIS  OF  THE  ELECTROMAGNETIC  FIELD  WITHIN  A  CIRCULAR 
WAVE  GUIDE  CONTAINING  AN  AXINLLY  NGN  UNI  FORM  PLASMA 

INI  RO DUCT  ION 

Tin*  purpose  of  this  chapter  is  to  relate  analytically  the  local  electromagnetic 
properties  of  the  plasma  with  measurement  s  of  certain  field  parameters  at  the 
boundary  of  the  wav.  guide  In  order  to  fix  ideas  consider  a  semi-inlinite  cir¬ 
cular  wave  guide  of  radius  R  whose  walls  constitute  a  perfect  conductoi.  Let 
this  guide  be  filled  with  a  plasma  having  properties  that  may  vary  arbitrarily 
only  in  the  axial  or  /.-direction.  The  configuration  so  described  is  a  simplified 
r  presentation  of  the  shock  tube  at  any  give!  time  during  the  propagation  of  the 
pressure  pulse  down  the  tube.  Ahead  of  the  disturbance  the  medium  is  yet 
unaffected,  and  alb  properties  are  uniform.  As  the  disturbance  is  trunsversed 
in  the  /.-direction,  the  pressure  and  temperature  of  the  medium  are  drastically 
altered,  as  ..re  its  properties.  Variations  of  properties  and  parameters  in  the 
radial  direction  are  neglected  here,  and  such  ctlects  are  discussed  elsewhere*. 

If  electromagnetic  radiation  is  introdu  ed  into  the  end  of  the  tube  towaid 
which  the  pressure  pulse  is  moving,  the  radiation  will  progress  down  the  tube 
through  multiple  reflections  from  the  perfectly  conducting  walls  of  the  tube 
until  it  encounters  the  ionized  medium  at  the  shock  front.  At  this  point  reflection 
and  absorption  of  the  radiation  will  begin  to  occur  and  this  phenomena  will  con¬ 
tinue  as  the  radiation  advances  into  the.  plasma  until  the  energy  is  completely 
returned  or  dissipated  Tne  resulting  field  at  any  point  within  the  tube  is  thus 

*  For  example,  see  General  Applied  Science  Laboratories,  Inc,.  Technical 
Report  No.  c55. 


a  composite  of  the  original  radiation  introduced  and  that  reflected  back,  from 
th  it  portion  of  the  plasma  extending  beyond  the  location  in  question.  Clearly 
the  field  pattern  will  be  characteristic  of  the  axial  distribution  of  electro¬ 
magnetic  properties  within  the  tube,  for  any  alteration  in  this  distribution 
would  result  in  different  degrees  of  absorption  and  reflection  at  each  point 
along  the  tube,  and  thus  the  resultant  field  would  be  changed. 

It  will  be  shown  in  this  chapter  that  for  the  simple  plasma  model  studied, 
namely  one  described  b\  the  two  independent  properties  of  collision  frequency 
v  and  plasma  frequency  the  two  characteristic  properties  may  be  measured 
indirectly  through  measurement  of  the  lime-average  ol  the  square  ol  the  radial 
component  of  the  electric  lie  Id  and  a  similar  average  for  the  tangential  component 
of  the  magnetic  field  and  the  axial  derivatives  of  these  quantities. 

In  view  of  the  fait  that  we  have  as:  umed  that  the  properties  are  uniform 
.it  any  given  axial  position,  such  measurements  may  be  carried  out  at  the  wall 
of  the  tube  where  r  1<  rather  than  by  the  use  of  probes  immersed  in  the  medium 
and  tending  to  alter  toe  environment.  In  order  to  measure  the  axial  variation, 
t  ie  detector  need  merely  traverse  the  tube  in  the  axial  direction,  or,  what  is 
equivalent  and  actually  used,  a  stationary  dete  tor  measures  the  desired  para¬ 
meters  as  the  pressure  pulse  moves  down  the  tube  carrying  with  it  the  charac¬ 
teristic  field  distribution. 

In  the  analysis  to  follow  the  rationalized  inks  system  of  units  is  used 
throughout,  and  the  nomenclature  is  the  conventional  one  and  is  tints  self- 


explanatory. 
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l.  FIELD  EQUATIONS 


Tlic  e lect  1  oma gnetii  Hold  within  the  wave  guide  described  in  the  preceding 


section  is  determined  from  the  Maxwell  equations 

DH 

dt 


V  x  E  -  -  p 


(D 


vx  n  -  j  1 1  iil 

pi 


(i> 


tlioir  associated  initial  conditions,  with  the  assumption  of  no  not  space  charge, 


V-  K  ..  0 

v  •  7i  o 


(3) 

(4) 


and  a  momentum  equation  relating  tn<  inrrent  density  to  the  electric  field: 


dt 


(5) 


The  plasma  is  nonmagnetic  and  unpolari/.ed  so  that  the  permeability  p  and 
permittivity  <  ar<  l  a  free  space  values.  1  ne  plasma  frequency  a>D  and  collision 
frequency  v  a  re  issumed  to  depend  upon  the  axial  coordinate  z  with  the  former 
defined  by  the  relation 


where  n  is  the  electron  number  density,  <>  its  charge,  and  m  its  mass. 

It  is  of  interest  to  note  that  for  a  transverse  electric  mode  (C7  -  0)  the 
current  given  by  (5)  is  also  transverse  and  is  thus  perpendicular  to  the  gradient 
of  the  electron  density  or  plasma  frequency.  This  behavior  is  consistent  with 
our  assumption  of  a  steady  distribution  of  i  lectron  density  depending  upon  the 
axiaL  coordinate  alone,  and  the  TE  mode  will  be  the  only  mode  considered  in 


this  analysis. 
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The  natural  system  of  coordinates  for  this  problem  is  the  circular 
cylindrical  system  (r,  t),  /.)  with  the  positive  /.-direction  being  taken  as 
opposite  to  the  direction  in  which  the  pressure  pulse  is  moving.  The  bound¬ 
ary  conditions  .ire  the  conventional  ones,  namely  that  the  tangential  compon¬ 
ent  of  the  electric  field  vanish  at  the  surface  of  the  perfect  conductor,  that 
the  field  he  finite  at  the  axis  of  the  tube,  and  that  tile  field  vanish  infinitely 
far  into  the  plasma.  Expressed  mathematically,  we  have, noting  that  we  have 
postal  ited  a  TE  mode  lor  which  Ez  =  d  everywhere, 


at 

r  =  R: 

K*  -  3 

(7) 

at 

r  =0; 

K  finite 

(8) 

as 

/.  — >  a>: 

E  — >  0;  f|  — 

0 

(9) 

It  may  be  shown  that  the  recessary  condition  that  the  normal  component  of  the 
magnetic  field  vanish  at  tin-  wall  is  automatic  ally  satisfied  by  t  ie  above  con¬ 
ditions. 

Following  the  usual  procedures,  the  pertinent  components  of  the  field 
satisfying  the  above  relutioi  s  may  be  shown  to  lie,  apart  from  a  factor  c110^, 

Er  F  Jn  ( V’>  (n  -  [)  (lJ) 

Ha  -  -i-  Kr  (vv'=  fj^i)  (U) 

«  |iu)  VV  UZ 

where  J,  is  the  Bessel  function  of  the  first  kind  of  order  n  and  the  constant 
Yj<_  may  have  any  of  the  discrete  values  given  by  the  roots  of 


The  function  vv( /,)  is  the  solution  of 


0 


(13) 


w 


00 


(13  a) 
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RELATING  AXIAL  VARIATION  OF  FIELD  TO 
AN”  (lld>AV 

It  may  be  shown  that  the  value  ol  ( 1”  )  ,  which  denotes 

A  V 

Ot  the  square  of  the  radial  component  of  tan  ,  lectric  field  is  Kivui  exactly  by 
the  relation 

(H,:) 


y  . ~~  the  time-average 


•W 


—  F.  E  * 
<  r  r 


(20) 


where  Lr  on  the  right  side  of  the  equation  i-  the  complex  expression  given 
by  equation  (Id)  or  (lb)  „,d  Er*  denotes  its  complex  conjugate  obtained  by 
replacing  i  with  -i  in  the  original  expression,  -imilails, 


(H9) 


w 


i  £ 

I  AA 


(21) 


Where  I Iy  on  the  right  sid.  of  the  equation  is  t  ie  complex  representation  of 
the  field  given  by  equation  (II)  or  (U).  in  o  ,  h  case  the  factor  e  iu,t  di  s 
and  f,irlhor’  be  can  .-a;  of  tie  manner  in  wl.i  O  tin  1 -dependence  e  nt 


tppea  rs 
ers,  namely 


in  the 

in*J 

lorm  < 

,  tins  factor  al 

so  di s 

app. 

<t  r 

Thus, 

t  ie  quantitie 

-  Kr  A* 

and 

11, 11 

1  th 

penc 

coordinates.  In 

p,. 

rtii  ul.tr, 

at 

r  =  R, 

Wl 

ha  vi* 

2  (  lA 

K  E* 

^  U“ 
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A  * 

1 

\V 

r  r 

r 

nk 

nk 

*i(]I  ) 

HI! 

A  A 

*  w 

\v’* 

W 

nk 

nk 

d  upon  only  the  radial  and  axial 
u.  ing  equations  (18)  mcl  (I)), 


„  * 


r  =  R 


(22) 


(23) 


The  importance  of  these  relations  lies  in  the  fact  that  the  axial  distribution 
of  the  quantities  ww+  and  w'w'*  may  thus  be*  determined  at  the  surface  of  the 
wave  guide  for  any  mode  through  measurement  "of  the 
and  (lly)  We  may  now  use  these  measurements  to  determine  K 


two  quantities  (  E 1 ) 

r  A  V 


1ZI 


3 .  DKTEKMINATION  OF  TI IE  SQUARE  OF  TIIK  LOCAL  COMPLEX^ 
REFRACTIVE  INI) EX,  ^ 

The  complex  retractive  index  *  is  related  to  the  function  w(  /-)  through 
the  differential  equation  (ll),  or  equivalently 


For  convenient  e, 


We  now  define  the  complex  function  F(/.)  by  the  left  .side  of 


(Z- 1),  namely 


F(  '■ )  = 


(25) 


We  desire  the 


real  and  imaginary  parts 


of  F ( / )  ,  and  thus  through  (Z5)  ,  deter¬ 


mine  k  . 

The  real  part  of  F(z).  using  (Z-i)  and  (lb)  i-  evidently  given  by 

w"  w1*  _vw\v*  w"* 

..Ref  “^r  ~~  w*  w\/f 

Noting  that  the  derivative  c.t  the  conjugate  is  equivalent  to  the  conjugate  of 


(Zb) 


the  derivative,  i  •  v.  , 

#  r 

w 1  s  w 

(Z7) 

and  using  the  identities 

(w  'w*)'  =  w"w*  -1  w'w* 

(Z8a) 

(w'Av )'  =  \V*"\V  t  w*  w' 

( Z  8b) 

Equation  (Zb)  may  be  written 

,  _  (ww*)"  -  Zw'w*' 

WW* 

and 


12Z 


F  =  w'w*  _  j_  (ww*) 
w\v*  Z  ww* 


(39) 


Referring  now  to  equations  ( Z Z)  and  (Zi),  we  have 

=  i  lRr),.W 

[:-'Kr  "  Kr^?  c:;)xv  ^  7^)w 


(30) 


and  thus  Fr  is  determined  in  terms  of  measured  quantities. 

Tiie  procedure  for  determining  the  imaginary  part  of  F  is  simiLer,  but 
not  as  straightforward.  Analogous  to  equation  (Z6),we  have 


ii  F.  = 

i 


w 

w 


It  tl 


V.'* 


l^w*  yv 

ww* 


(31) 


and  by  an  obvious  transformation 
U*  I  in  u 


Fi 


(33) 


Tiie  function  w'/w  may  be  determined  in  the  following  way.  First  we 
note  that 

W'  4  *'*  .  (WW*)'  _ 


W  U*  ww* 

where  a  is  a  measurable  function.  Also,  another  measurable  function  is 
w'w* ' 


(33) 


ww* 


s  b 


(3d, 


Solving  (33)  and  (3!)  simultaneously  for  w'/w,  we  ha 


r 

i  -4  —  =  _L  i  t 

_  - 

w*  a 

1  -  \ 

^  (w'w* ')  (ww*r 

w  Z 

i 

k  a2  Z  (ww*) 

(  W  \\  *  )  'j  ” 

(35) 


1Z3 


Returning  to  equation  (3Z)  we  now  have,  substituting  from  equation  (35), 


Im  7i_  =  1  i  (ww*)  ,  (Ww*  )  (ww*)  _  ( 

w  Z  (  ww*)  Y  £(ww*)'^ 


Fi  =  1 


(ww*) 


(w'w*1)  (ww*) 

4  — k - - —  i 


*(  w  w*  p^j  - 


Z  ww* 


(36) 


(37) 


Using  again  equations  (ZZ)  and  (Z3),  we  have  linally 


F.  =  ' 
i 


•^"b  <"«  )rW  <Kr  Kv 


.v. 


Z  (":r'Av 


(38) 


The  choice  of  sign  in  cquatioi  (38),  as  will  be  seen  in  the  next  section, 
is  made  to  yield  a.  positive  collision  frequency. 

The  real  and  imaginary  parts  of  «  or  k  “  ma>  now  be  determined  directly 
from  equation  (Z5) , 


k  -  =  11  F  • 


-7  'k  5  K  +  1 


(39) 


T  hu  s , 


Rex" 


,  .L 

Im  k  =  T  ~  F. 

w“  1 


(40) 


(41) 
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4.  PLASMA  FREQUENCY  AND  COLLISION  FREQUENCY 

The  desired  quantities  u>p  and  v,  or  what  is  equivalent,  <J>  and 
readily  found.  From  the  definition  of  k1  ,  equation  (15),  we  have 


4>*  \  _  j 

1+  y  /  1 1  o  " 


Util 


Solving  for  y  and  cp  ,  we  have 


I  _ 

R-  t 


0"  =  (  1  -  R)  +  — -  -  (l  R)  (l  t  y 2 ) 

(1-1  ) 


Usim;  equations  (10)  md  (41),  we  have 


>  0 


and 


-i 

Y^)  (1  ’  /) 


The  plasma  properties  are  thus  determined  hv  the  measured 


41  are  now 

(42) 

(4  3) 

(44) 

(45) 

(16) 

quai  lilies 


th rough  the  use  of  equations  (30)  and  (3K)  for  Fr  and  Fj  . 
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